Pattern formation in 2-frequency forced parametric waves 
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We present an experimental investigation of superlattice patterns generated on the surface of a 
fluid via parametric forcing with 2 commensurate frequencies. The spatio-temporal behavior of 4 
qualitatively different types of superlattice patterns is described in detail. These states are generated 
via a number of different 3-wave resonant interactions. They occur either as symmetry-breaking 
bifurcations of hexagonal patterns composed of a single unstable mode or via nonlinear interactions 
between the two primary unstable modes generated by the two forcing frequencies. A coherent 
picture of these states together with the phase space in which they appear is presented. In addition, 
we describe a number of new superlattice states generated by 4-wave interactions that arise when 
symmetry constraints rule out 3-wave resonances. 

PACS numbers: 47.54. +r, 47.35. +i, 47.20Gv 



I. BACKGROUND 

Patterns are ubiquitous in the world around us. The 
word "pattern" describes an order, regularity, or a simple 
mathematical description, which can be found either in a 
natural or man-made system. Patterns often result from 
self-interactions of driven nonlinear systems. Naturally, 
the first patterns to be scientifically analyzed were the 
simplest ones, which can be described by few mathemat- 
ical variables. However, in recent years we have learned 
to recognize and categorize patterns in systems that were 
assumed to be formless, devoid of any order. Perhaps 
the most obvious characteristics of these systems are the 
multiple length and time scales that can be present si- 
multaneously. One of the most important mechanisms 
to explain such phenomena is the nonlinear resonant in- 
teraction between the different modes that are excited in 
these systems. In these interactions two or more waves 
can interact to form "new" waves. These waves have a 
wavelength and frequency which is the sum or difference 
of the basic waves. The system's energy can then be 
transferred between these modes or dissipated at differ- 
ent scales. The purpose of the work described in this 
paper was to explore this paradigm in the experimental 
study of a simple controlled system: the parametric ex- 
citation of waves on the surface of a fluid (the Faraday 
system). 

The general form of the external acceleration applied 
to the system is given by: 

g(t) = A[cos(x) cos(muj t) + sin(x) cos(nw ^ + 4 1 )] (1) 

This spatially uniform vertical excitation preserves the 
system's spatial symmetries, while modifying its tempo- 
ral ones. 

As first noted by Faraday, sinusoidal acceleration (in 
the direction of gravity) of a fluid layer with angular fre- 
quency lu induces a pattern, having a wavenumber k(oj), 
on the fluid surface. Whereas waves excited by a single- 
frequency have been studied extensively over the last 4 



decades, the response of the system for multi-frequency 
excitations has only recently begun to be investigated. 
Single-frequency driving can produce patterns of differ- 
ent symmetries. Patterns consisting of rolls, squares, 
hexagons and 8,10,12-fold quasi-patterns have been ex- 
perimentally observed in |j]-H]. The secondary instabili- 
ties of these different patterns involve complicated states 
that display either spatiotemporal chaos (| J8|j|,g| , trans- 
verse amplitude modulations |l(J] or various defects pJlCJ] 
that break the patterns' initial global symmetry. The use 
of multiple-frequency driving enables us to study the in- 
teractions of different excited modes in a controlled way, 
as each excitation frequency can linearly excite a well- 
defined wavenumber. In this way, we hope to be able 
to slowly unfold the system's underlying behavior and, 
thereby unravel the fundamental mechanisms that de- 
scribe the waves' interactions. 

The purpose of this paper is to provide a coherent 
overview of the wide variety of nonlinear states that re- 
sult from 2-frequency forcing. We will provide detailed 
descriptions of the spatial and temporal behavior of these 
states. In doing so, we will provide a characterization 
of these nonlinear states - depicting both the resonant 
mechanisms and symmetry constraints giving rise to their 
formation. 

This paper is organized in the following fashion. In 
section Q we will briefly describe the theoretical and ex- 
perimental work that has, to date, been performed in this 
system. The experimental apparatus and measurement 
techniques used in our measurements will then be de- 
scribed in section D. We will then present, in section III, 
an overview of the phase diagram together with a brief 
description of the different types of superlattice states 
observed. Each type of superlattice, together with the 
mechanisms that form it, will then be described in detail 
in the subsequent sections. A codimension two point, 
at Xci exists in this system where both externally driven 
modes simultaneously become linearly unstable. Subhar- 
monic Superlattice States (SSS), which bifurcate from an 



initial hexagonal state far from Xc will be described in 
section IV . We will then progress to the region of phase 
space in the vicinity of Xc- Three different types of super- 
lattice patterns will be described in sections |V|, [VJ, and 
VII. All of the above superlattice patterns result from 



different types of 3-wave resonant interactions. We will 



conclude with section VIII in which a number of super- 



lattice states generated by 4-wave resonant interactions 
are described. We will show that these states can occur 
when 3-wave interactions arc forbidden. 



A. Notation 

The notation conventions used throughout this paper 
are as follows. The driving function is specified in Eq. (Q). 
To avoid confusion, we will always specify the driving fre- 
quency ratio used in Eq. (gj) by the ratio m/n where m 
and n are the two co-prime integers that describe the 
two frequencies u)% — muo and 101 = hloq. We will always 
assume that n > m. In our notation, k\ and ki refer 
to the wavenumbers excited, respectively, by the driving 
frequencies u>\ and ui^. The angle x m Eq. ([j]) describes 
the relative mixing between the two modes and the an- 
gle <f> describes their phase difference, where the relevant 
range is < <f> < 2irm/n . cos(x) and sin(x) are some- 
times 1 11 12] replaced by the mixing coefficients r and 
1-r. 

In the following sections we will frequently characterize 
eigenmodes by their temporal parity. Since the parity of 
m and n is important we will use the notation odd/even 
(even/odd) to described the classes of driving where m 
(n) is odd and n (m) is even. Odd/odd describes driving 
where both m and n are odd. A state whose temporal 
response has a fundamental frequency of luq/2 will be 
denoted as a "subharmonic" state, whereas a "harmonic" 
state is one with a fundamental frequency of luq. 

We shall use the following units. The total ampli- 
tude A, appearing in Eq. ([l]), is measured in units of 
g = 981cm 2 / sec. The fluid's kinematic viscosity, v is 
measured in cS (0.01cm 2 /s), and the depth of the fluid 
layer h, is measured in cm. In many cases, we will iden- 
tify, for simplicity, angular frequencies (e.g. to) with the 
corresponding temporal ones (e.g. u)/(2tt)). Where nec- 
essary, angular or temporal frequencies will be explicitly 
denoted. Additional notation will be defined as needed. 



B. Linear analysis 

The linear stability analysis of the problem was per- 
formed by Besson et al. [n3[ by numerically solving the 
linearized Navier-Stokes equation via an extension of the 
technique developed by Kumar and Tuckerman Mm for 
single frequency excitations. As in the single frequency 
case, the acceleration-wavenumber plane is characterized 



by alternating tongues corresponding to the acceleration 
at which a given wavenumber becomes linearly unstable. 
For mixing angles x — 0° and x = 90° the tongue struc- 
ture of single frequency forcing with either mujQ or tiluq 
driving is obtained. Increasing x from results in the ap- 
pearance of additional tongues whose dominant frequen- 
cies are spaced uio/2 apart, since the system's basic fre- 
quency is then ujq. Each odd numbered tongue possesses 
a subharmonic temporal dependence composed of only 
frequencies u>o(p + 1/2) whereas even- numbered tongues 
are temporally harmonic i.e. composed of frequencies 
pu)o (where p is any whole number) . Although the time 
dependence of each tongue is given by an infinite series, 
the dominant frequency of thep th odd (even) tongue cor- 
responds to (p+l/2)u>o (pwo). Generally, within the crit- 
ical tongues the dominant frequencies nuJo/2 or muo/2 
have an order of magnitude greater amplitude than the 
other components. Thus, the temporal response of the 
critical modes for m/n =odd/odd driving ratios is always 
subharmonic. For odd/even (even/odd) the k\ "tongue" 
has a subharmonic (harmonic) response while the &2 has 
the opposite parity. 

The system's critical wavenumber, k c , corresponds to 
the lowest acceleration, a c , at which the flat, feature- 
less state loses stability. At a critical value of x = Xc 
a codimension 2 point exists where two tongues having 
wavenumbers k\ and ki corresponding, respectively, to 
mwo/2 and nwo/2, simultaneously become unstable. For 
X far from Xc the critical wavenumbers are close to the 
values of k c obtained for single-frequency experiments. 
Near the codimension 2 point k c tends to differ from the 
corresponding single- frequency value by up to 10%. The 
calculated values of both a c and k c are in excellent agree- 
ment with experiments 13 la]. While modes other than 
the critical ones are linearly damped, we will see that 
they can play an important role in nonlinear wave inter- 
actions. 



C. Experiments with two- frequency forcing 

In the case of single-frequency driving, the subhar- 
monic time dependence prohibits quadratic terms in the 
"amplitude" equations describing the nonlinear interac- 
tions between the amplitudes of the excited modes. How- 
ever, when using two-frequency driving with odd/even or 
even/odd parity, both harmonic and subharmonic tem- 
poral responses are possible. When one of the driving 
frequency components is dominant, one can consider the 
smaller component as a perturbation that breaks the sys- 
tem's temporal subharmonic symmetry. The reflection 
invariance of the corresponding set of coupled ampli- 
tude equations is then broken and, generically, quadratic 
terms can appear. These quadratic terms are important 
since they enable three-wave interactions between differ- 
ent modes. 



Edwards and Fauve were the first to study the two- 
frequency driven Faraday instability Jl6|— |l8[] . They chose 
to focus most of their study on 4:5 driving although 
they also explored other ratios (such as 6/7, 4/7, 8/9 
and 3/5). These experiments used a relatively viscous 
fluid and small fluid layer height in order to minimize 
lateral boundary effects. As the viscosity of the fluid was 
rather high (y = 100 cS) stripe patterns occurred for 
single-frequency driving. The phase space as presented 
in Jig ] for (even/odd) 4:5 driving can be divided into two 
parts: The harmonic (subharmonic) part where the k\ 
(fc 2 ) wavevector is dominant and the leading temporal 
term has the frequency of Aloq/2 (5wo/2). In the har- 
monic region, in place of the stripe patterns of wavenum- 
ber k\ appearing for pure 4wo driving, a first order transi- 
tion to hexagons occurs for \ > 10° • In the subharmonic 
region, striped patterns with wave number fc 2 are ob- 
served until the near vicinity of the codimension 2 point 
at Xc- The first-order transition to the hexagonal state 
results from the quadratic interactions mentioned above. 

In the neighborhood of the codimension 2 point, a tem- 
porally harmonic, twelvefold symmetric quasi-periodic 
pattern was observed. These states appeared for only a 
small range of 4> {<j> ~ 75° ±5°). They evolved, via a first- 
order bifurcation, from either the flat zero-amplitude 
state or the subharmonic striped patterns. In the hys- 
teretic region of these states, Edwards and Fauve also 
observed solitary axial waves that originated from the 
quasi-periodic pattern. Arbell and Fineberg |19] have 
shown that these highly locallized waves are related to 
"oscillons" states observed [E0( in driven granular sys- 
tems. 

Muller 1T3] later conducted two-frequency forcing ex- 
periments using a 1/2 driving ratio. These experiments 
were performed near Xc in shallow fluid layers in various 
regions of the x ~~ 4> phase space. Both temporally sub- 
harmonic hexagon and triangle patterns were observed. 
Triangular patterns are formed when the spatial phase 
associated with each of the excited eigenmodes differs 
from 0° or 180°. Muller showed that amplitude equations 
with both cubic and quintic terms (applying to tempo- 
rally subharmonic waves), can form triangular patterns. 
In contrast, amplitude equations with quadratic terms 
have only stable hexagonal solutions. These experiments 
were later modelled by Zhang and Vinals p/M using a 
quasipotential approach. Experimentally, Muller showed 
that the addition of a third small-amplitude forcing fre- 
quency (which is equivalent to perturbatively breaking 
the system's parity) could stabilize either the hexagonal 
or triangular states. 

More recent experimental studies in 2 frequency-forced 
systems were performed by two groups, Kudrolli, Pcir 
and Gollub J2l) and Arbell and Fineberg 0@|| . 
These studies were conducted both in the near vicinity 
and far from Xc- They revealed a number of qualitatively 
new, superlattice-type states in which new scales, not di- 



rectly introduced via the external forcing, were evident. 

In regions of phase space that can be relatively far 
from Xc superlattice states were observed as secondary 
bifurcations from the harmonic (rowo/2) hexagonal states 
that occur for odd/even or even/odd driving ratios. The 
primary hexagonal symmetry with wavenumber k c is 
broken by additional modes with wavenumbers q < k c 
whose temporal response possesses an mwo/4 component. 
These states include "SL-II" states observed for 4:5 driv- 
ing by Kudrolli et al. plf and the Subharmonic Superlat- 
tice ( "SSS" ) states observed for a large number of driving 
ratios by Arbell et al. E5 |. 

A second type of superlattice state is observed in the 
near vicinity of the codimension 2 point, once again on 
the side dominated by the harmonic driving component 
when even/odd or odd/even forcing is used. Two variants 
of these states coined "SL-I" [glj and double hexagonal 
states ("DHS") |li| have been observed. These states 
can be described by the superposition of two hexagonal 
sets of wavevectors of magnitude k c . The two sets of six 
wavevectors are oriented at an angle 9 r ~ 22° to each 
other. This specific angle is not arbitrarily chosen. The 
sum and difference vectors between the two wavevector 
sets produce a sublattice, spanned by the smaller differ- 
ence wavevectors. When the two sets of wavevectors are 
oriented at specific angles of 6? r , the sublattice formed 
by the difference vectors becomes commensurate with 
the two hexagonal lattices. This structure is one of the 
generic possibilities that were anticipated on the basis 
of symmetry arguments proposed by Silber and Proctor 

A qualitatively different type of superlattice occurs in 
the vicinity of Xc for all driving ratio parities. These 
states, coined 2-mode superlattices (2MS) p5i, are the 
most general of the superlattice states observed. They 
are formed by an interaction of the two linearly excited 
modes (k\ and fc 2 ) with a third linearly-damped slaved 
mode that is nonlinearly excited. The angle between k\ 
and ki is chosen by the following resonance mechanism: 
the vector difference fc 2 — k\ produces a third wavevector 
k^. The magnitude of k% is determined by the dispersion 
relation, w(fc), for the difference frequency W3 = w 2 — u)%. 

An additional type of superlattice state has been ob- 
served for 2/3 and 4/5 driving in the vicinity of Xc p2[ . 
This state, which appears in place of the 2MS state, con- 
sists of a rhomboid pattern that is formed by the simple 
nonlinear resonance: fc 2 — fc' 2 = k\ where |fc 2 | = |fc' 2 |. 
When the coupling angle 6 (9 = cos _1 [fc 2 • fc 2 /fc|]) be- 
tween fc 2 and k'2 is tuned to a value of 9 ~ 27r/n, 2-n 
fold quasicrystalline patterns are naturally formed. 



D. Model equations and nonlinear analysis 

Generally, two methods have been used to study the 
Faraday instability with two-frequency driving. The first 
uses simple model systems that yield qualitative insights 
regarding the behavior of the Faraday system. These use 
general assumptions based mainly on symmetry consid- 
erations. The second method is to start from the full 
nonlinear set of equations that describe the system, em- 
ploy carefully chosen approximations, and derive a set 
of equations that describes the behavior of the system 
based on the real physical parameters. Both methods 
have yielded valuable insights. 



1. Model equations 

The observation of quasi-crystalline patterns ( "quasi- 
patterns" ) generated using two-frequency driving by Ed- 
wards and Fauve |lq| and via single frequency driving 
by Binks et al. |24) provided a motivation to find model 
equations that display similar behavior. Muller [g5| first 
considered a system of N coupled Landau equations with 
cubic nonlinear terms. These equations could be writ- 
ten as the gradient of a Lyaponov functional. Muller 
showed, by minimization of this functional, that regular 
N-fold patterns of different symmetries can be stable. 
Pattern selection depended on the value of the nonlinear 
coefficients coupling the linearly degenerate modes. This 
mechanism may be related to both the appearance of 
quasi-patterns in single frequency Faraday experiments 
and to the quasi-patterns observed in the harmonic re- 
gion of two-frequency Faraday systems with even/odd 
driving. 

A second mechanism that can create quasi-patterns 
is related to quadratic interactions between degenerate 
nonlinear modes. Muller proposed that a quadratic non- 
linearity, generating the triad interaction: k 2 — k 2 = k\ 
flfel = |Aj |) could also lead to quasi-crystalline pat- 
terns. The angle between k 2 and k 2 is tunable by the 
ratio ki/k 2 , with resonant angles 9 = 45°, 9 = 36°, 
9 = 30° for 8-fold, 10-fold and 12-fold quasi-patterns 
respectively. These states were observed in a system 
of two coupled Swift-Hohenberg equations, each with a 
different unstable wavenumber. Frisch and Sonnino ]2q| 
also observed subcritical 10-fold symmetrical patterns in 
coupled Swift-Hohenberg equations. This state was nu- 
merically shown to be stable even when the dynamics 
are not derived from a free energy functional. In ad- 
dition to states with N-fold symmetry, nonsymmetric, 
rhomboidal patterns were also seen to be stable for some 
parameter values. Later, both the rhomboidal patterns 
and resonant quasipatterns resulting from the above in- 
teractions were observed experimentally by Arbell and 
Fineberg p|]. 



Lifshitz and Petrich |2j modeled the two-frequency 
Faraday system with a single generalized Swift- 
Hohenberg type equation for a single real field u(x,y). 
This model is simpler than the coupled equations used 
by Muller, Frisch and Sonino, and Newell and Pomeau 
p8[ . The model equation used was rotationally invariant 
with two built-in critical wavenumbers. The equation 
contained a quadratic term that both broke the system's 
up-down symmetry and allowed triad wave interactions. 
Stable striped, hexagonal and 12-fold symmetric patterns 
were observed for different values of the control param- 
eter. In addition to these N-fold symmetric states, a 
compressed hexagon state, similar to the rhombic/stripe 
pattern described by Muller, was observed. 

The above model systems suggest that the existence of 
two unstable wavenumbers together with the possibility 
of triad interactions (provided by quadratic terms) is a 
sufficient condition for the formation of quasi-periodic 
patterns. Another common feature of these models is 
the existence of distinct regions of phase space in which 
patterns that lack N-fold symmetry are stable. 



2. Nonlinear analysis 

In contrast to the simplified model systems described 
above, Zhang and Vinals [O derived a description of 
the system's dynamics from the governing equations for 
the 2-frequency Faraday problem. To this end, they ap- 
plied the quasi-potential approach developed for single- 
frequency study p9| to the problem of two-frequency 
driving. This approach is strictly valid in the limits of 
weak dissipation and infinite fluid depth. 

To compare their results with Muller's experimental 
results, Zhang et al. analyzed the special case of 1/2 
driving in depth. They first used the linearized equa- 
tion to study the location of the codimension two point, 
Xc as a function of the phase difference 4>. The results 
were in qualitative agreement with the experiments. The 
discrepancies were attributed to the high damping used 
in Muller's experiment, which was outside the region of 
validity of the theory. It is interesting to note that the 
dependence of \ c on is a special feature of 1/2 driving 
and does not occur for other driving combinations. 

Zhang and Vinals then, using a multiple scales ap- 
proach, derived standing wave amplitude equations. This 
weakly nonlinear analysis assumed that the system was 
far from the codimension 2 point, so that a single tem- 
poral mode dominated the dynamics. For the case of 1/2 
driving, they first obtained a prediction for the relative 
magnitudes of the different Fourier components of the 
weakly nonlinear temporal response of the fluid surface. 
Then, assuming N degenerate modes, the coupled ampli- 
tude equations describing these modes were derived. In 
contrast to |25| . the function f3(9ij) coupling the i th and 
jth m0( i eSj was computed from the physical parameters 



of the problem. One interesting result of this calculation 
was that the phase difference (j> can have a strong effect 
on the coupling function [3 - and thereby a strong effect 
on the nonlinear pattern selection. The relative stability 
of different N-fold nonlinear states was then calculated 
by minimization of a Lyapunov functional, as in pq ] and 
p7j . Semi-quantitative agreement with the regions of 
stability observed experimentally in [Q for different pat- 
terns (squares, hexagons/triangles and quasi-patterns of 
various orders) observed in the subharmonic region of the 
4> — x phase space was obtained. 

An important result of this work was that it suggested 
a new type of physical mechanism that governs the se- 
lection process. These calculations indicated that the 
value of the coupling coefficients was strongly influenced 
by triad interactions between the linearly excited modes 
(corresponding to the dominant excitation frequency) 
and the linearly stable modes corresponding to the second 
excitation frequency. Resonant coupling to these latter 
modes served, in the region of phase space far from \c to 
enhance the effective damping - as energy transferred to 
these "slaved" modes is more efficiently dissipated. Far 
from x c , states that can not couple to the slaved modes 
are then preferred by the system. As we shall later see, 
however, in the vicinity of \c resonant triad coupling to 
the linearly stable, slaved modes provides one of the main 
mechanisms for the rich variety of nonlinear states ob- 
served. 

Silber and Skeldon f3G] were the first to theoretically 
study the two-frequency Faraday system in the vicinity 
of the codimension 2 point. This study pointed out the 
importance of accounting for the temporal symmetries of 
the system. Silber and Skeldon focused on forcing ratios, 
m/n, having either odd/even or even/odd parities, where 
interactions between harmonic and subharmonic waves 
may occur. 

As shown in |[Ll| , resonant mode interactions greatly 
affect the mode coupling function j3(0). Using normal 
form analysis, Silber and Skeldon showed that triad res- 
onances k\ ± k[ — ki (where k\ — k[) are only pos- 
sible when the temporal mode corresponding to fe is 
harmonic. When ki has a subharmonic temporal de- 
pendence, quadratic terms in normal form equations can 
be eliminated (3^3^] - thereby decoupling the harmonic 
modes from the subharmonic ones. This can be simply 
understood since the product of two linear eigenfunc- 
tions (resulting from a quadratic interaction term) re- 
sults in the addition of their temporal phases. The sum 
of two harmonic or subharmonic temporal phases can- 
not produce a subharmonic one, therefore two modes of 
like parity cannot couple quadratically to a subharmonic 
state. Silber and Skeldon [M went on to demonstrate the 
above, by calculating the amplitude equations for both 
1/2 (odd/even) or 2/3 (even/odd) driving by means of 
the quasi-potential approximation used in pdj| . 



The existence or suppression of three wave resonances 
can have a significant effect on the qualitative features of 
the phase diagram. When one is far from \c we have seen 
pd| that three-wave resonant coupling influences pattern 
selection by enhancing dissipation via the coupling to a 
heavily damped (slaved) mode. In this case, resonant tri- 
ads are strongly suppressed. On the other hand, when in 
the near vicinity of Xc resonant coupling to nearly un- 
stable linear modes can occur. Moreover, if these modes 
undergo a first order bifurcation, their growth will be 
(to first order) unchecked and resonant coupling to them 
may have a very significant effect on the spatial-temporal 
behavior of the system. 

Silber, Topaz and Skeldon (g3[ have recently demon- 
strated the importance of resonant coupling to slaved 
modes for 6/7 forcing near \c- Using the quasi-potential 
approximation flllfl they showed how weakly damped lin- 
ear modes with wavenumbers K < k c , quadratically cou- 
ple to the unstable modes to create the "SL-F states 
observed by Kudrolli et al. 21 1. In this case, the critical 
wavevectors \ki\ = k c could be constructed from a com- 
mensurate hexagonal sublattice of wavevectors, Ki, such 
that ki = qK\ + pK2- The "SL-I" states are a particular 
case where (p, q) = (±2, ±3) with k c /K — \fl. This par- 
ticular coupling was made possible by the existence of a 
weakly damped linear "tongue" with a wavenumber close 
to K. The "SL-I" state was not observed for 2/3 forcing 
since, for this forcing ratio, no linear tongues near this 
resonance exist since there are no additional harmonic 
modes with K < k c (Ki must be harmonic by []30f| as 
they result from the vector difference of two ki modes) . 

Recent work by Tze, Rucklidge, Hoyle and Silber 
|34f has shown that the "SL-II" states observed by Ku- 
drolli et al. [gl| may be understood as resulting from 
a symmetry-breaking bifurcation of an initial hexagonal 
symmetry. Study of the possible invariant subgroups of 
the original Dq + Z 2 symmetry characterizing hexago- 
nal standing waves revealed a number of possible solu- 
tion branches. One of these corresponds to the spatial 
symmetry of temporally-averaged "SL-II" type states. 
Depending on the normal form coefficients, 5 additional 
possible solution branches were predicted. It remains to 
be seen whether these other branches are experimentally 
observed. 



II. EXPERIMENTAL SYSTEM 

Our experimental system consisted of a shallow fluid 
layer, laterally bounded by a plastic sidewall and 
mounted on a computer-controlled mechanical shaker. A 
1 cm thick, black-anodized aluminum plate of 14.4 cm 
diameter supported the fluid from below. This plate was 
machined to a 10/im flatness. The mechanical shaker 
used (either Unholtz— Dickie model 5PM or VTS model 



100) provided vertical accelerations ranging from to 
15g. The cell acceleration, regulated to within O.Olg, 
was monitored continuously by a calibrated accelerome- 
ter (Silicon Designs, INC 1210L-010) attached directly 
to the armature of the shaker. A feedback mechanism 
was used to control and stabilize the amplitude A, mix- 
ing angle x, and phase </> to desired values. 

Most of our experiments were conducted with 
Dow Corning 200 silicone oil of different viscosities 
(DC200/10, DC200/20, DC200/50 and DC200/100). Sil- 
icone oil has a typical density of 0.95 gr/cm 3 and surface 
tension of 21.5 dyne/cm. This fluid is Newtonian for 
the viscosity range of 1-100 cS. Since the fluid viscos- 
ity is highly temperature dependent, stabilizing the fluid 
temperature was important. A stable fluid temperature 
of 30 ± 0.05°C was used in all experiments. Resultant 
viscosity variations were less than 0.04 cS. A number of 
experiments were also performed using TKO-77 vacuum 
pump fluid with viscosities ranging between 221 cS at 33° 
and 184 cS At 30°. Both Dow Corning 200 and TKO-77 
have very low vapor pressures so there was no need to 
seal the cell against evaporation. The results of our ex- 
periments showed no dependence on the particular type 
of fluid used. 

Our experiments were performed at frequencies be- 
tween 20-150 Hz. The selection of the frequency was 
influenced by the aspect ratio of the patterns and the 
shaker's maximum acceleration and stroke. Frequency 
selection was also influenced, to a lesser extent, by limita- 
tions of the imaging and laser probe technique. Typically 
the aspect ratio between the cell diameter, L, and the 
wavelength A, was between 5 < L/X < 50. The maximal 
driving frequency of 150 Hz was governed by a c , which 
increases with increasing u>. The shaker's maximal stroke 
(2.5 cm peak to peak) and boundary mode quantization 
at small aspect ratios dictated the lower frequency limit. 



A. Boundary Conditions 

The lateral boundary conditions of the experimental 
cell can have an important effect on the waves excited by 
the system. In our experiments we attempted to mini- 
mize the role of the sidewalls. A circular shape for the 
lateral boundary was chosen. This ensured that no par- 
ticular pattern was preferred. This is especially signif- 
icant when the system is only slightly dissipative (e.g 
low viscosity fluids and/or large fluid depth) . For a more 
highly dissipative system (e.g high viscosity fluids and/or 
shallow fluids), the boundary's shape does not influence 
the symmetry of the excited pattern JOI . 



i 



FIG. 1. (top) A profile of the container used is shown: a 
Delrin circular boundary (gray) is attached to the bottom 
plate (black). The boundary consists of a vertical section of 
height ho and an inclined section at an angle of a = 20° 
chosen to allow the surface of the liquid used, Dow Corning 
200, to be at zero contact angle with the rim (liquid in light 
gray). Rings with ho of 1 mm, 1.5 mm, 2.5 mm, 4 mm and 
5 mm were used, h could be changed continuously by adding 
small amounts of fluid with a calibrated pipette. 
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FIG. 2. (a) A schematic diagram of the experimental sys- 
tem. The imaging system consisted of a cylindrical screen, 
concentric with the experimental cell, which was illuminated 
by halogen lamps arranged in a circle. A CCD camera was 
mounted on the cylinder axis above the fluid surface, (b) The 
cylinder's illumination intensity was varied as a function of 
the height above the fluid. At each point on the fluid surface 
the local slope reflects only a single point from the cylinder 
surface into the CCD. Since the lighting provides a unique 
intensity at each height along the cylinder, the intensity re- 
flected by each point is uniquely mapped to the projection of 
the fluid surface's slope in the direction of the cylinder axis. 

As discussed by Douady (1990) pj, an additional ef- 
fect of sidewalls is the possible emission of waves (menis- 
cus waves) from the lateral boundaries. These waves are 
forced at the driving frequency via forced height varia- 
tions of the meniscus formed at the contact line between 
the fluid and lateral boundaries. Meniscus waves have no 
threshold and can therefore mask the instability thresh- 
old of parametrically forced waves. To minimize this ef- 
fect, our system's lateral boundaries (shown schemati- 
cally in Fig pi) were sloped at an angle conjugate to the 
contact angle between the fluid and the material (Delrin) 
from which the lateral rings were constructed. (A slope 



of a = 20° was used for the Dow Corning 200 silicone oil.) 
In this way, we ensured that the static fluid interface was 
nearly flat. 

An additional advantage of sloped lateral boundaries is 
the elimination of reflected waves by impedance match- 
ing. Since the instability threshold increases with de- 
creasing h, a gradual decrease (sloping sides) in fluid 
depth increases the effective local threshold at the larger 
radii to far beyond a c . Since the typical height of the 
fluid layer in the sloped region was only 0.1 — 0.8mm, 
parametric waves could not be excited and any meniscus 
waves emanating from the wall were strongly damped. 
In practice, this boundary condition combined with the 
fluid viscosities and depths used enabled us to obtain val- 
ues of a c within 2% of the calculated values |^,[l3|,|36| for 
a system of infinite lateral extent. 



B. Visualization 

1. Imaging from above 

To visualize the fluid surface, we employed a novel 
imaging technique. The imaging system is schematically 
shown in Figure |2jb. The experimental cell was illumi- 
nated by a tall cylindrical screen whose axis was concen- 
tric with center of the cell. The screen was illuminated 
from below by a ring of 12 small lamps. As a result, the 
light intensity along the screen varied as a function of the 
height above the fluid. A CCD camera was mounted on 
the cylinder axis, 1.4m above the cell. At each point on 
the fluid surface (see Figure ||b), the local slope reflects 
only a single point from the screen onto the CCD. Since 
the lighting provides a unique intensity at each height 
along the cylinder, the intensity reflected by each point 
on the fluid surface is uniquely mapped to the projection 
of its slope on the cylinder axis. We used the CCD's high 
speed shuttering mode (1/1000 sec) to obtain instanta- 
neous images of the fluid surface. 

Two methods of triggering were used to control the 
CCD camera. The first method employed a trigger signal 
that was synchronized with the driving. This signal both 
reset the camera and initiated acquisition of the video 
frame at a desired phase relative to the driving signal. To 
observe slow changes in the patterns over long times, slow 
trigger rates that were commensurate with the driving 
frequency were used. The short-term behavior of a state 
in its different temporal phases was studied by the use 
of slightly incommensurate trigger rates. This allowed 
nearly continuous acquisition of the different temporal 
phases of a given state without the need for very high- 
speed acquisition. 

Our imaging technique, although providing quantita- 
tive information, does not directly yield the surface wave 
height function h(x,y). The imaging yields a grayscale 



image, I(x,y), that is approximately the absolute value 
of the gradient of the height function i.e. 

I(x, y) = J(d x h(x,y)) 2 + (d y h(x,y))* (2) 

One must then work backwards from I(x, y) to deter- 
mine the function h(x,y). This is done by inputting an 
assumed state into Eq. H and comparing the computed 
pattern to the state observed. By iteration it is possi- 
ble to arrive at fairly good estimates of h(x,y). Some 
examples are presented in Figure |3|. 
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FIG. 3. Using the method described in the text we calcu- 
lated model images (upper) of a simple hexagonal state in 
its peak state (a) and crater state (b) and of a square state 
(c). The corresponding experimental patterns are shown (bot- 
tom). To model the surface waves we assumed an asymmetry 
between up hexagons and down hexagons due to the the fact 
that the hexagon's amplitude was large compared to the small 
layer height. 

For high-amplitude states, visualization from the side 
was sometimes used. This was performed by illuminating 
from the side and placing a video camera in the horizontal 
plane of the plate at the height of the system's lateral 
boundary. This configuration enabled direct quantitative 
measurements of the wave amplitudes adjacent to the 
lateral boundaries of the cell. An additional advantage 
of this imaging was that it allowed us to simultaneously 
view both the lower plate's vertical movement together 
with the wave's motion. In this way their relative phase 
could be directly measured. 



C. High resolution temporal measurements 

Most previous studies of the Faraday system have used 
imaging techniques that mainly yield information about 
the symmetries of the states observed. Time domain 
information was generally extracted by the use of stro- 
boscopic lighting at the frequency of the shaker. This 
technique measures the lowest frequency in which a state 
oscillates but yields no information about higher harmon- 
ics. Two experimental studies used a laser beam probe to 
extract information about the surface waves. Douady O 



used a laser beam deflected at an angle from the surface 
waves and reflected on a screen to study the amplitude of 
the waves. Gollub and Simonelli |J7| used a laser beam 
deflected from the surface waves onto a photo detector to 
obtain temporal information. This technique, however, 
could not give precise information about the slope of the 
fluid surface since an average intensity reflected from a 
single "cell" of the pattern was measured. 

We combined these two methods by imaging the re- 
flection at the fluid surface of a highly focused laser onto 
a position sensitive detector (UDT SL20). This method 
yielded an accuracy of 1 — 5% in the surface slope and 
a temporal resolution of 0.02 mS. Good signal to noise 
was obtained by amplitude-modulating the laser signal 
and then deconvolving the resultant signal of the posi- 
tion sensitive detector (PSD). The temporal response of 
the system was only limited by the maximal sampling 
frequency of the PSD voltage. 



III. OVERVIEW OF THE PHASE DIAGRAM 

The phase space of the two-frequency Faraday system 
is very large. Besides physical parameters such as fluid 
layer height and viscosity, one also has to set the driving 
parameters. Introducing two driving components with 
different frequencies, different amplitudes and a nontriv- 
ial relative phase makes the task of choosing a working 
regime and the relevant dimcnsionless combinations of 
parameters a difficult one. 

We have chosen to focus on the simplest commensurate 
driving ratios as a first step. The m/n ratio combinations 
used in our experiments were numerous: 1/2, 2/3, 2/3, 
2/5, 2/7, 3/4, 3/5, 3/7, 4/5, 4/7, 5/6, 5/7, 5/8, 6/7, 6/11, 
41/60, 40/59, 21/50. Most of our detailed experiments 
were performed with no phase difference between the two 
frequency components (0 = in Eq. (|l|)). For system pa- 
rameters that were seen to excite special patterns, scans 
of <j) were made. Changing (f> was found to be crucial for 
the existence of some of the states and of no relevance to 
others. Phase diagrames were constructed by fixing the 
mixing angle x an d increasing the amplitude A until a 
state of droplet ejection was reached. 

Two typical two-frequency phase diagrams for 
even/odd driving are presented in Figure 0. In sin- 
gle frequency experiments rolls, squares, hexagons and 
quasipatterns of different symmetries are known to exist 
depending on the viscosity, height, frequency and am- 
plitude above the threshold. In our parameter regime, 
the dominant structure in regimes dominated by a single 
frequency is squares for low viscosities (y — 8.7 — 23 
cS) and both squares and rolls for higher viscosities 
{v = 47 — 87 cS). As in other studies of two-frequency 
driving |13,nq,[l7|, two main regions, dominated by either 
k\ or fc2, exist. Each of these regions has not only a dif- 
ferent wavelength but also different temporal behavior. 



The fci and k^ dominated regions have a strong response 
at uji/2 and W2/2 respectively. This is true for all values 
of m/n. The temporal response, however, consists of ad- 
ditional frequency components that depend on the ratio 
m/n. When both m and n are odd (odd/odd driving) 
the excited surface modes in both regions of phase space 
have only subharmonic components (i.e. (j>+1/2)ujq with 
p an integer) . In the case of odd/even (even/odd) driv- 
ing, the k\ dominated region is temporally subharmonic 
(harmonic) while the k^ dominated regime is temporally 
harmonic (subharmonic). At the critical value \ — Xci 
a codimension 2 point exists where both wavenumbers 
are simultaneously linearly unstable. The k\ dominated 
region occurs for x < Xc while the ki dominated region 
occurs for x > Xc- The interaction of k\ and k^ leads to a 
variety of different nonlinear states in the vicinity of Xc- 
Before describing these states, we will first describe the 
effects of two-frequency driving in the two main regions 
far from Xc- 

In the regions dominated by harmonic states we have 
found a number of non-trivial states that bifurcate from 
single-mode hexagonal states while breaking both their 
spatial and temporal symmetries. These symmetry- 
breaking bifurcations can even occur when x is quite 
small. The subharmonic superlattice state (SSS) (Figure 
(bottom)) is an example of such a symmetry breaking 
state. SSS states are formed in the k\ dominated regime 
with Xc ~ 5° > X > 10°, when the primary hexagonal 
state's symmetry with harmonic temporal behavior is 
broken by an additional set of wavevectors of magnitude 
q < k\ with subharmonic temporal behavior (frequency 
wi/4) with respect to the primary lu±/2 frequency. We 
have observed two main types of SSS states, which dif- 
fer from each other in the orientation, magnitude and 
number of q wavevectors. For example, the SSS state 
shown in Figure |], (SSS type I) breaks the initial hexag- 
onal symmetry by the introduction of two wavevectors 
q, which are parallel to two of the three initial wavevec- 
tors k\. The magnitude of q in this case is fei/2 which 
yields a simple resonance condition q + q = k±. SSS have 
been observed only for even/odd driving for all of the 
m/n combinations listed above. The different SSS types 
and the mechanisms that form them will be described in 



detail in Sec. IV 



For even/odd driving, the effect of two frequency driv- 
ing on the pattern formation in the ki dominated region 
is quite different than in the k\ dominated region. In 
the ki region, square symmetry dominates at threshold 
from x — 90° to x ~ Xc- Only in the vicinity of Xc do 
we see the effects of the two-frequency driving on the 
patterns formed. It is interesting to note that although 
theory predicts that hexagons are preferred for harmonic 
response and squares for subharmonic response (see Sec. 
PJ), we have observed square symmetry in large parts of 
the harmonic region for odd/even driving (x > Xc)- 




Unlocked Hexagonal 2MS 2kR DHS 

FIG. 4. Two typical phase diagrams of two-frequency ex- 
periments obtained for 2:3 driving with different system pa- 
rameters: loo — 20 Hz and h = 0.155 cm (left) and too — 25 
Hz and h — 0.2 cm (right). In both experiments v — 23 
cS and <f> — 0° . Square regions exist in the near vicinity of 
single frequency forcing. The square symmetry dominates in 
the subharmonic regime to near \ a while throughout most 
of the harmonic region the hexagonal symmetry dominates. 
In the vicinity of the co-dimension 2 point we observe three 
new states that exist for many combinations of the driving 
ratio. These states are 2 mode superlattices (2MS), with un- 
derlying square or hexagonal symmetries and spatially and 
temporally unlocked states ("Unlocked" states). A resonant 
state that consists of a rhomboid unit cell (2kR) was also ob- 
served (right). Unlike the 2MS and unlocked states, which 
appear for many different driving ratios (odd/odd, odd/even 
and even/odd), this state was observed for only for 2/3 and 
4/5 driving. In the harmonic region of phase space where 
hexagons are initially dominant, a second bifurcation occurs 
to either temporally subharmonic states (Subharmonic Su- 
perlattice State, SSS) or high amplitude waves (as well as, 
at times, localized "oscillon" waves) that appear on a double 
hexagonal superlattice (DHS) . Symbols in the phase diagram 
describe measured transitions for fixed \- (bottom) Typical 
photographs of these states. 

Let us now briefly describe the patterns formed for 
even/odd driving in the vicinity of \ c . Starting with the 
k\ dominated (\ < Xc) region, two types of patterns are 
observed near \ c . One pattern, which we call the "Dou- 
ble Hexagonal State" (DHS) is formed by two sets of 
hexagonally arranged wavevectors (of length fei) with a 
finite angle, a, between them. In the phase space shown 
in Figure [|, a ~ 22°. In contrast to the SSS, this state 
does not break the temporal symmetry of the harmonic 
hexagon state. Depending on various system parame- 
ters, DHS are sometimes formed by a first-order bifur- 
cation. Perhaps their most outstanding characteristic is 
their very high amplitude. The surface wave maxima 
can reach amplitudes much higher than the fluid layer's 



height. In JL9| we have shown how the DHS states can 
form oscillons, a highly localized large-amplitude nonlin- 
ear state that has been observed ]l9| , pO|J3"g| ] in a number 
of periodically driven systems. 

A special case of the DHS occurs for a = 30° , whence 
one obtains 12-fold quasipatterns such as first observed 
in Mm. We have seen the formation of such patterns 
for 4/5 driving in the same region where the DHS with 
a = 22° appears for 2/3 driving. 

Let us now move to the vicinity of Xc both on the bor- 
der of the k\ region and within the ki region. Here two 
linearly unstable eigenvectors with different magnitudes 
can be concurrently excited. There are numerous possi- 
ble configurations in which such a system can organize 
itself. Four different classes of mixed mode states were 
found to exist. These will be described in detail in Sec- 
tion 



VI 



• Two mode superlattices (2MS) These states formed 
by the interaction of the dominant mode (e.g. fo) 
with its original symmetry and the weaker mode 
(e.g. k\), which breaks the symmetry of the origi- 
nal pattern. The symmetry of the dominant mode 
can be either square or hexagonal depending on 
the proximity of the nearest primary state in phase 
space. In Figure we present two types of 2MS 
modes, a square 2MS obtained for x > Xc and a 
hexagonal 2MS obtained for x < Xc- The temporal 
behavior of the 2MS contains both the wi/2 and 
W2/2 frequency responses and always includes the 
subharmonic frequency of ljq/2. These states ap- 
pear for all types of driving (odd/even, even/odd 
and odd/odd) although the precise structure of 
phase space depends on the driving ratio used. 

• Unlocked States Between the square and hexago- 
nal 2MS states an intermediate region exists where 
both k\ and ki appear but no well-defined sym- 
metry or spatial mode locking is observed. Thus, 
no long range correlations in either space or time 
exist. The basic time scale of the surface waves is 
T — Att/lvo but the pattern changes its local struc- 
ture over time scales of order (10 2 — 10 3 )T. 

• Rhomboidal states Changing h, v or luq can lead 
to qualitative changes in the phase space. Figure 
shows two different phase diagrams obtained for 
2/3 driving. The only difference between the two 
diagrams is the fluid layer height and the value of 
basic frequency luq. Lioubashevski et al. [B9[ de- 
scribe how the dimensionless number S/h, defined 
by the ratio between the effective boundary layer 
depth, 5 = yJv/ojQ, and the fluid height, h, af- 
fects single mode states selected by the system. 
For a certain range of S/h rhomboidal patterns re- 
place the 2MS and unlocked states (see Figure 
bottom). These states couple two wavevectors of 



length &2 with one wavevector of length ki . These 
wavevectors evolve spontaneously from two circles 
of linearly degenerate states. We will show that for 
special parameter values, n-fold quasipatterns can 
naturally evolve from the rhomboidal structures. 

In addition to the states described above, in Section 



VIII we will show examples of a number of other reso- 
nant structures that are formed for different system pa- 
rameters. The richness of this system allows one to ob- 
serve a wide variety of different resonant patterns. We 
will show some common traits of these resonant selection 
mechanisms that can lead to a more comprehensive un- 
derstanding of resonant interactions in pattern forming 
systems. 



IV. SUBHARMONIC SUPERLATTICE STATES 

(SSS) 
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FIG. 5. For 2/3 driving we observe two types of temporally 
subharmonic superlattices, the SSS-I (a) SSS-II (b). Experi- 
mental images (left) together with their spatial power spectra 
(middle) and simulated images (right) are shown. In (a) we 
see that a small vector of magnitude q = K/2 (q]\K) breaks 
the hexagonal symmetry formed by primary wavenumbers Ki, 
where Ki = k c . In (b) the primary hexagonal symmetry is 
broken by wavevectors q = K/y/3 located at a 30° angle rel- 
ative to the primary wavevectors. Both patterns can exist in 
regions far from \ c , where the higher odd-frequency compo- 
nent is weak (see Figure ^) . The patterns obey the resonance 
conditions 2q=k c (SSS-I) and q + q> = K (SSS-II). 

Hexagonal patterns can be formed in two-frequency 
experiments when the basic subharmonic temporal in- 
variance is broken and quadratic nonlinear terms appear 
in the amplitude equations describing the system. These 
terms can occur in regions of the two-frequency phase 
space that have a harmonic response. By increasing the 
amplitude of the driving, the primary hexagonal symme- 
try can be broken and new stable structures appear even 
in regions where the second externally driven mode can- 
not be excited at all. In this section we will describe the 
characteristics of these symmetry breaking patterns. 



In the notation used in this section, K is the linearly 
unstable wavevector, excited io\ = mwo, that charac- 
terizes the primary pattern, q is an additional smaller 
wavevector that appears in states that bifurcate from the 
primary pattern. Since in each pattern there is a degen- 
eracy in the direction of the wavevectors, we will use an 
index to number the different wavevectors of the same 
magnitude e.g. Ki i = 1, 2 — 

Sub-harmonic superlattice states occur over a wide 
range of x m t he two-frequency phase diagram where the 
lower frequency, oji, is dominant. SSS were observed for 
most even/odd frequency ratios tested, but were not seen 
for odd/even or odd/odd ratios. Since the first SSS type 
state was found |15| additional variants have been ob- 
served for different experimental parameters. Both types 
of SSS pattern, SSS-I (Figure |a) and SSS-II (Figure |b) 
are shown with their corresponding spatial power spec- 
tra. Both states share the following characteristics: 

• The basic temporal dependence of both SSS types 
is subharmonic with respect to the the primary in- 
stability i.e. the system has a strong response at 
muJo/A. 

• These states are secondary bifurcations of tempo- 
rally harmonic hexagonal states. 

• Both states occur in the same basic region of phase 
space (for < x < Xc)- 

• In both states &2 wavevectors are not observed. 

The different states, in general, occur for different fluid 
parameters such as fluid viscosity, fluid depth and uji/lu2 
ratios. The most obvious differences can be seen in Fig- 
ure |5[ Examining their spatial power spectra, we see that 
while both SSS-I and SSS-II's spectra are constructed by 
three evenly spaced wavevector pairs, Ki. The spectra in- 
clude additional sets of wavevectors, ql (where \q%\ < K), 
of smaller magnitude. The orientations and magnitudes 
of ql differ in SSS-I and SSS-II states. 

In SSS-I ql are oriented solely along (some or all of the) 
axes defined by the Ki. The lengths of the qi in the SSS-I 
are usually K/2, but in some cases the qi ^ K/2. 

SSS-II states are formed by a set of qi which always 
consists of three evenly spaced wavevector pairs of magni- 
tude q — K/^/3 arranged at an angle of 30° relative to the 
direction of the Ki triad. The magnitude and orientation 
of qi yield the simple resonance condition ql + qi = Ki . 
While the SSS-I have been observed for a wide range of 
driving ratios (2/3, 4/7, 2/5, 2/7, 4/5), the SSS-II have 
only been seen for 2/3 and 4/5 driving. Both types of SSS 
states have recently been identified as representations of 
different invariant subgroups when hexagonal symmetry 
is broken |0|Q. 

In the next subsections we will present a detailed de- 
scription of each of the SSS types together with a mech- 
anism that can explain their formation. 
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A. Subharmonic Superlattice type I (SSS— I) 

The spatio-temporal behavior of the SSST state can 
be modelled by a simple equation for the surface height 
function: 

1 3 

h(r, t) = cos(-mwot) 2_\ ^ cos (^i ' ^+ a i) + 

i=\ 

1 M 

cos(-muj t + 7) y^ Bj cosjcjl ■ r + fa) (3) 



2=1 



with 



q-WKt, A =0,90° 



(4) 



where M is the number of axes with broken symmetry, 
7 is the temporal phase difference between the two sets 
of modes and oti and Pi are the respective spatial phase 
of the Ki and ql components. Because the pattern has 
hexagonal symmetry we assume on = 0° and not on = 
2ir/3 as in patterns having triangular symmetry [p.2| . Eq. 
summarizes the most important features of the SSST 
in a compact way. We will now present experimental 
evidence for the validity of this equation and describe 
the relevance of each term (Ai, Bi, 7, 0i and M) in the 
experimentally observed states. 
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FIG. 6. Typical temporal sequences of an SSS-I state (a) 
taken at constant values of the driving parameters for the fre- 
quency ratio 40/60 Hz. The spatial Fourier spectra are com- 
posed solely of wavenumbers K, corresponding to Wi/2, and 
q — K/2. The locations of the peaks are fixed in the different 
temporal phases. For the system parameters ojo/(27t) — 20 
Hz, zx = 23 cS, 4> — 0°, and h = 0.155 cm the hexagonal sym- 
metry is broken in two directions by the K/2 vectors. Circles 
of radii K and K/2 are drawn in (b). The relative intensities 
of the different wavevectors can be seen; in (a) qi are nearly 
absent while in (c) their intensities are almost equal to the 
Ki . The symmetry breaking is also revealed in the intensities 
of the primary hexagonal vectors, as can be seen in (c), where 
the two strong intensity wavevectors are enclosed in a square 
and the weak one in diamonds. 






FIG. 7. Images (top) and power spectra (bottom) of an 
SSS-I type state with broken symmetry in one, (a) two, (b) 
and three (c) directions. The circles indicate the primary 
hexagon wavevector's magnitude, K. The symmetry break- 
ing vectors of magnitude q = K/2 are enclosed in squares. 
All images were obtained for the same system parameters of 
v = 23 cS (40/70 Hz and h = 0.155 cm, x = 43° and <p = 0°). 
This pattern is not stable and drifts slowly (order of 10-60 
seconds) between these three states. 

In Figure || we present a sequence of SSS-I states taken 
at different times for constant values of the driving pa- 
rameters. Although the state's appearance changes with 
time, their spatial Fourier spectra reveal that the state 
results from the interaction of two specific spatial scales; 
the primary wavenumber \K\ = \ki\ that is excited by the 
u>i frequency component, and its spatial sub-harmonic, 
q = K/2. The SSS spectra show that while the Ki have 
6-fold symmetry, two q\ with relatively large amplitude 
and a third smaller amplitude qi have broken this sym- 
metry. The amplitudes of the qi vary with time. Within 
the temporal phases shown in Figure oa the qi amplitudes 
have little power while in Figure oc the qi are stronger 
than the Ki components. This behavior is reflected by 
the temporal phase, 7, in Eq. a. The symmetry break- 
ing is reflected by the relative strengths of both the qi 
and Ki. It is clear that the two K wavevectors enclosed 
in squares, have different strength than the wavevector 
enclosed in diamonds. This symmetry-breaking is also 
seen in the relative power of the qi wavevectors. 

As demonstrated in Figure 0, SSS states can have bro- 
ken symmetry in one, two or three directions (M — 1, 2, 3 
in Eq. [3|). In most cases, a specific number of symme- 
try breaking directions was selected. However, for the 
experiment described in Figure |7], the pattern drifted 
slowly (on the scale of seconds) between the three possi- 
ble symmetry-broken states. 

The transition between the primary hexagonal pattern 
and the SSS states is (within ss 1%) non-hysteretic and 
occurs via a circular front that propagates slowly inward 
from the plate's lateral boundaries. The process is con- 
tinuous and reversible (see Figure ph. Before continuing 
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the study of the spatial characteristics of the SSS-I pat- 
tern let us digress briefly and examine the temporal be- 
havior of the system using the laser probe method (see 
Sec. |ilC|). 



bilize certain patterns such as the 2kR and oscillon states 

MM. 




FIG. 8. The transition between Hexagonal and the SSS 
state is displayed for a typical experiment (40/60). This is 
a gradual process in which the basic symmetry of a perfect 
hexagonal pattern (a) is broken first at the circular boundary 
of the cell. As the amplitude of the external forcing increases, 
the area of broken symmetry grows inwards (b,c) until the 
SSS state fills the entire cell (d). This process can also occur 
in the reverse direction. 

Typical time series of the x and y components of the 
fluid surface gradient at a single point are presented in 
Figure g for three different accelerations. These describe 
the temporal behavior of hexagons at threshold (a), de- 
veloped hexagons (b) and SSS-I (c). At threshold, the 
response is harmonic with respect to the total period of 
2-k/uq. The response is strongest at the frequency of 
uji/2. Increasing the driving amplitude results in a bifur- 
cation to a state with a strong superharmonic response 
at the frequency of iO\ . This phenomenon also occurs for 
single frequency experiments. 

A further increase of the driving amplitude yields a 
second bifurcation. In this bifurcation the temporal re- 
sponse becomes subharmonic with respect to the period 
of 2tt/ojq. As can be seen in Figure |9| the superharmonic 
component does not disappear and can be quite strong. 
It is important to note that 2/n driving has some spe- 
cial relations between the various frequencies that are not 
present for higher order driving (such as 4/5, 4/7. ..6/7...). 
For 2/n driving Wi/2 = wo- The temporal response in 
the harmonic region has the same frequency as the com- 
mon frequency lu . In other driving ratios, such as 4/5 
driving, the subharmonic is 4cjo/2 = 2ujq where the com- 
mon frequency is u>q. It is possible that such a relation 
can enhance certain resonant mechanisms and help sta- 
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FIG. 9. The time dependence of the SSS-I state as stud- 
ied by the reflection of a laser by the surface waves. In this 
experiment, for system parameters wo/(27r) = 22 Hz, v — 23 
cS, h — 0.2 cm, <f> = 0° and \ — 36.3°, increasing the driving 
amplitude from 2.5 g (a) to 2.9 g (b) and 3.7 g (c) results in 
(a) Low amplitude hexagons (b) developed hexagons and (c) 
SSS-I states. In (a) both d x h and d y h have the same peaks 
in their power spectrum. Increased driving amplitude yields 
a different temporal response in the two directions x and y. 
This phenomenon also occurs in single frequency experiments 
where squares are dominant and may be a general feature of 
the Faraday instability in viscous fluids. In the SSS-I state 
(c) a subharmonic temporal response at uq/2 occurs. 

Wagner et al. Q studied a two mode system gener- 
ated by single-frequency excitation at a bicritical point 
where both harmonic and subharmonic tongues become 
unstable. They describe two different superlattice states 
that exist in the transition region between subharmonic 
squares to harmonic hexagonal states. In the first super- 
lattice state, the square symmetry is broken by a small 
wavevector kr> that is equal to kn — ks, where ku, ks are, 
respectively, the primary hexagonal and square wavevec- 
tors. The relative phase of the symmetry breaking kr> 
mode compared to the primary square ks mode can ei- 
ther be 0° or 90° according to the sign of the nonlin- 
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ear coefficient in the amplitude equation for the symme- 
try breaking mode. The experimental observation shows 
that the phase selected for the first transition state is 90°. 
A second superlattice state J42| was observed that is sim- 
ilar to the SSS-II states found in our experiments. The 
spatial phase difference, /3, between the harmonic hexago- 
nal mode and the subharmonic symmetry-breaking mode 
was found to be 0° although from amplitude equation 
considerations 90° is also a possible solution. 
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FIG. 10. A simple one-dimensional superposition of 
a harmonic function (bottom) cos(27rx/A) with its sub- 
harmonic can have two basic combinations. The "dis- 
placive" mode shown results from the superposition of 
cos(27ra;/A) + bcos(-Kx/\ + n/2) (top) (b = 2.2 is arbitrar- 
ily chosen), whereas a "modulational" mode results from the 
superposition of cos(27ra;/ A) + &cos(7ra;/A) (middle). It can be 
seen that in the "displacive" mode the distance between local 
maxima is either d or 2A — d with A — d <C A, while in the 
"modulational" mode the distance between the local maxima 
remains the same as in the primary mode, A. 

We will now describe in detail the effects of this phase 
difference on SSS-I states with broken symmetry. Let us 
first consider the effect of adding a K/2 modulation in 
one dimension. In Figure |l0| we can see that adding the 
K/2 mode with /% = 0° results in a modulation of the 
original cosine form. The waves' local maxima remain at 
the same spatial location but their amplitudes are now 
modulated with a An/K periodicity resulting in one large 
peak followed by a smaller peak. In contrast to this, 
adding the K/2 mode with j3i = 90° results in a different 
effect. While all of the maxima have the same ampli- 
tude, their spatial locations are modulated with a 4n/K 
periodicity. This effect |42| is called "displacive". The 
spatial effect in images can be quite pronounced. This is 
demonstrated in Figure O where we present a compar- 
ison between the "displacive" (left) and "modulational" 
(right) effects with spatially subharmonic patterns and 
symmetry breaking along one direction (top) or three di- 
rections (bottom) . The effects of symmetry breaking in 
three directions are more complicated but the qualitative 
effect remains. We find that SSS-I states are always of 
the displacive nature ((3 — 90°). This can be seen by 
comparing the experimental images shown in Figure M 
with the simulated displacive patterns shown in Figure 

U 



It is interesting to note that the "modulational" pat- 
tern shown in Figure O (bottom right) was seen in an ex- 
perimental study of the single-frequency Faraday system 
in a viscoelastic fluid Gj3] . When changing the driving fre- 
quency, Wagner et al. found both a harmonic region for 
low driving frequencies and a subharmonic response for 
higher ones. In the vicinity of the transition frequency, 
Wagner observed a hexagonal superlattice composed of 
both the subharmonic and harmonic wavevectors. Since 
the subharmonic wavevector is exactly half of the har- 
monic one, the simple resonance of k$ + k$ = kn is re- 
tained (where the S index stands for subharmonic and the 
H for harmonic) . Although the power spectra of this su- 
perhexagon described by Wagner are similar to the SSS-I 
power spectra (in the case of three direction symmetry 
breaking), the real space patterns are different, indicat- 
ing a 0° spatial phase difference between k$ and kn, in 
contrast to the /3j = 90° seen for SSS-I. 
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FIG. 11. "Displacive" (left) and "modulational" patterns 
(right) are shown for both symmetry-breaking in 1 (top) and 
3 directions (bottom). The pattern shown is simulated using 
our imaging model applied to Eq. (W) with M=l (M=3) for 1 
direction (3 directions) and with /3; = -k/2 ( /?, = 0) for the 
"displacive" ("modulational") patterns. 2:3 forcing was used 
with amplitudes of all modes taken to be equal (Ai — Bi). 
All SSS-I patterns were found to be displacive in character 
(compare to Figure m). 

The combination of the spatial displacive mechanism 
and the subharmonic temporal dependence of the K/2 
mode results in an interesting "jittering" effect in time. 
If we consider Eq. || we see that when t — > t + 2-k/loq the 
first part of l.h.s of Eq. is invariant whereas the second 
part changes sign due to the different time dependence. 
In our model calculation the temporal displacement of 
2tt /uiq is equivalent to a spatial displacement of 2ix /K in 
each of the symmetry-breaking directions. Images pho- 
tographed at time intervals of 27t/cj appear to jitter be- 
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tween these two displaced images, as demonstrated in and 



Figure 12 
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FIG. 12. The "jittering" effect in 2d. Experimental (top) 
and model (bottom) images of the SSS-I states at two tempo- 
ral phases separated by a 2tt/ioo interval. The white crossed 
lines are fixed in both frames. The two pictures demonstrate 
that the time displacement is equivalent to a spatial displace- 
ment of 2tt/K in the symmetry-breaking directions. 



B. Subharmonic Superlattice type II (SSS— II) 

As mentioned above, a qualitatively different type of 
SSS type pattern, SSS-II, has been observed. Increasing 
both viscosity and height but keeping the dimensionless 
parameter S/h constant (by changing u>q), results in a 
different symmetry breaking scenario. Though the SSS- 
II appears in the same region of phase space as the SSS-I 
and shares its subharmonic temporal behavior, it is qual- 
itatively different from the SSS-I. The excited vectors ql 
are now aligned at an angle of tt/6 relative to the vectors 
Ki (in contrast to qi\\Ki in type SSS-I). The symmetry 
breaking wavevectors are of magnitude \qi\ = \Ki\/\f?> 
and these states exhibit no spatial symmetry breaking. 
For each Ki there is a corresponding ql. 

Our experiments suggest that SSS-II can be described 
by: 

1 3 

h(r, t) = cos(-mLOot) 2_, ^i cos(Ki ■ r + oti) + 

i=l 

1 3 

cos(-mwoi) 2_\ ^i cos (<7i ■ r + A) (5) 



<7i =g(-l,0);g"2 



,1 V3 

9(2'— i;* 



A -y/3. K 



Looking at this state at different temporal phases, the 
effect of the qi can be easily seen at some phases (Figure 
|l3|a) while, at others, only the hexagonal symmetry is 
apparent (Figure 13c). This symmetry is a representa- 
tion of one possible invariant subgroup when hexagonal 
symmetries are broken [EOpl^Mj] . 
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FIG. 13. Images (a,b left) with corresponding power spec- 
tra (a,b right) of different temporal phases of an SSS-II state 
observed for 2/3 driving in the harmonic region of phase space 
at v = 47 cS (iuo/(2n) = 25 Hz and h — 0.2 cm). In contrast 
to the SSS-I pattern, where the primary hexagonal symme- 
try is broken by wave vectors parallel to Ki, SSS-II feature 
a different symmetry breaking. A second hexagonal lattice of 
smaller magnitude wave vectors forms, obeying the resonance 
condition q{ + q% — K . At certain temporal phases the ef- 
fect can be easily seen (a) while at others only the hexagonal 
symmetry is dominant (b). 

In the 1/2 driving experiments performed by Mullcr 
p2| a transition between hexagonal and triangle patterns 
was observed that corresponds to a change of J^. oti in 
Eq. (S) from the value of to 37r/2 (where the second 
term of the R.H.S. of Eq. (pi) is zero). A similar phe- 
nomenon can take place with SSS-II states. As shown 
in Figure [lj, a variant of the SSS-II states occurs with 
a rotational symmetry of 27r/3 in contrast to the 27r/6 
that is typical for hexagon patterns. These SSS-II vari- 
ants are observed for 2/3 driving when relatively higher 
frequencies (u>o > iOHz) are used. Using our imaging 
model, we find that these patterns are formed by the 
same resonance as the SSS-II but when spatial phase 
angles of the <fi modes are fy = ir/2 for i — 1,2,3 or 
^2i&i = 37r/2. These SSS-II states have symmetries 
that are similar to those of the time-averaged symme- 
tries and wavenumbers of the "SL2" states observed in 
plf . The instantaneous images of SL2 states, however, 
are more reminiscent of SSS-I states (as shown in Fig- 
ure 0). Patterns similar to this SSS-II variant have re- 
cently been observed in a forced ferrofluid system Eq] 
in the vicinity of a bicritical point where harmonic and 
subharmonic solutions collide. Muller shows that when 
considering a harmonic region amplitude equation with 
quadratic terms for a single wavenumber model, only so- 
lutions where J^ /3j = are possible. It appears that 
the second set of equations for the q wavenumber modes, 
which are temporally subharmonic, does allow the exis- 



14 



tence of modes with J^i Pi — 37r/2 solutions. The mecha- 
nism that selects the spatial phase in multi-mode systems 
still must be clarified. 
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FIG. 14. For 2/3 driving and v = 0.47cS, a SSS-II type 
pattern with a different spatial dependence occurs. Two 
phases of this state are shown in (a,b left) and a simulation 
of these pattern are shown in (a,b right) . These patterns are 
formed when the spatial phases of the smaller wavevectors, 
(<f;), are shifted by n/2 relative to Ki. These states were 
observed for higher frequencies (ujq > 35 Hz) than those for 
which SSS-II having a zero phase shift are observed. 

SSS-II type patterns have also been observed in the ex- 
perimental study of optical pattern formation in sodium 
vapor |46[ . In the next section we will address the ques- 
tion of a selection criterion between the SSS-I and SSS-II 

states. 



C. Selection between SSS-I and SSS-II 

Both SSS-I and SSS-II break the temporal symmetry 
of the initial hexagon state by temporal period-doubling 
to a basic frequency of mwo/4. The two states, however, 
differ in their spatial behavior. For the case of SSS-I 
states, the excited wavevectors, ql, are both parallel to 
the linearly unstable wavevectors Kj, and, in many cases 
q = K/2. In contrast, the excited wavevectors in SSS-II 
states are rotated by tt/6 relative to Ki. What mecha- 
nism governs the selection of both the two different states 
and the values of q that are excited? 

As Silber, Topaz and Skeldon have suggested [p3[ , 
the symmetry-breaking wavevectors of harmonic patterns 
may correspond to minima of linearly stable tongues that 
can be excited via nonlinear coupling to the Ki modes. 
Since the subharmonic frequency, mwo/4, is excited by all 
of the SSS states, the linearly stable tongue with a dom- 
inant muJo/4: frequency would be a likely candidate to 
be selected. The wavenumber, q, corresponding to these 
waves, can be well-approximated by the linear dispersion 
relation q = k(muio/A). Note that q is not constant for a 
given value of u)q, but can be strongly dependent on the 
parameters v and h. 

Let us now examine the following premise. The system 
will generically prefer to undergo spatial period doubling 
to q = K/2. If, however, q(muJo/A) is close to a wavevec- 
tor q with a magnitude that is substantially different than 
K/2, one possible solution of the system is to lock to ei- 
ther SSS-I or SSS-II patterns with q ^ K/2. If q is near 
K/\/3, the system will lock to this value, thereby fulfill- 



ing the spatial resonance condition: ql + qj = Ki. SSS-II 
patterns will then occur. Values of q sufficiently far from 
either spatial resonance will result in an SSS-I state with 
q 7^ K/2. Such a state is presented in Figure [15] where 
symmetry breaking occurs in either one direction (Figure 
|l5|a,b ) or three (Figure hoc) but the symmetry-breaking 
wavevector's magnitude was not K/2. Instead, vectors 
parallel to Ki, with magnitudes q =/= K/2 and K — q are 
observed, with, empirically, q ~ 0.6k c . 




FIG. 15. Images (top) and corresponding power spectra 
(bottom) of SSS-I superlattices observed for 2/3 driving in 
which q 7^ K/2. Circles of a radii K (outer) and q, corre- 
sponding to the muo/4 frequency (inner) are drawn. The 
primary hexagonal symmetry is broken either in one direc- 
tion (a,b) or in three directions (c). (a,b) show the same 
state at different temporal phases. The relative intensities of 
q and K — qvaxy in the different phases. All experiments were 
with v — 23 cS. (a,b) for h — 0.2 cm and (c) for h — 0.25 
cm. The basic frequency was u>o/(2tt) = 35 Hz for (a,b) and 
wo/(27r) = 30 Hz for (c). The value of q/K is « 0.6 for all 
experiments. 

This premise is checked in Figure uM, where we plot 
the value of the ratio between the experimentally mea- 
sured values of q and the wavenumber computed for sin- 
gle frequency driving using mwo/4 with K for param- 
eters where different SSS-type patterns were observed. 
The plot shows that for both SSS-II and the SSS-I with 
q =/= K/2 the approximation q — fc(mwo/4) is correct 
to within 4%. In SSS-I states where q = K/2 is the 
symmetry breaking mode, the ratio q/k(mcoo/A) varies 
systematically between 77% — 85%. This suggests that 
the q = K/2 resonance is strong enough to induce this 
"locking" or detuning of q. 

For 2/3 driving we observed the appearance of SSS-I 
states at lower liquid layer depth and the SSS-II at higher 
depths. As h is reduced, we found that instead of an 
abrupt transition between the SSS-I and SSS-II states, 
both types of symmetry breaking can occur simultane- 
ously. As shown in Figure O, the primary hexagonal 
symmetry of this state is broken by two 6-fold sets of 
wavevectors. 

One set corresponds to SSS-I with magnitude qo = 
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K/2 while the other set corresponds to SSS-II, with q\ = 
K/y/3. In real space the characteristic pattern of SSS-II 
(Figure |l7p) is broken by superimposed stripes in one 
direction. This state has both modulative and displacive 
effects implying that the two symmetry-breaking modes 
retain their respective spatial phase characteristics. 
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FIG. 16. The measured wavenumber, q m eas, of the sym- 
metry breaking wavevectors in SSS-type patterns, normal- 
ized by the linear wavenumber calculated for the frequency 
muio/A as a function of the measured critical wavenumber of 
the primary hexagon pattern, K. The symbols correspond to 
different types of SSS patterns; SSS-I patterns with q = K/2 
(circles), SSS-I with q ^ K/2 (diamonds) and SSS-II patterns 
with q — K/\fj, (squares). The data indicate that the first 
type of pattern (q = K/2) is generally preferred by the sys- 
tem unless the value of q(muio/4) is either close to q = K/\/2> 
or sufficiently far from either of these preferred modes. 

As demonstrated in Figure [Hj the mechanism that 
forms both the SSS-I and the SSS-II patterns depends 
on a "slaved" linear eigenmode. In the 2/3 driving de- 
scribed above the only existent linearly stable subhar- 
monic tongue occurs for loq/2 — mtvo/^- In contrast 
to 2/3 driving, even/odd driving ratios with m/n val- 
ues such as 4/5, 6/7 etc. possess linearly stable tongues 
whose dominant frequencies differ from mwo/4. For ex- 
ample, in 60/75 Hz driving, besides the 15 Hz (= mtoo/4) 
response observed for the SSS state, a response at 7.5 
Hz (= ujq/2) is also possible. Recent theoretical work 
by Silber et al. |33[ has suggested that these additional 
slaved modes can influence the character of the selected 
nonlinear state. We find that these additional slaved 
modes can indeed appear. In 4/5 driving typical SSS- 
I states with a frequency response of mwo/4 are observed 
far from \c- Near Xc, both the 12-fold temporally har- 
monic quasicrystalline states first described by Edwards 
and Fauve Jig ] , as well as a state composed of a cascade of 
symmetry-breaking bifurcations occur. This second type 
of pattern, which has a subharmonic time dependence, is 
shown in Figure [lqa. 




FIG. 17. Reducing h from h = 0.33 cm where only SSS-II 
was observed to h = 0.31 cm results in a state where the two 
types of SSS patterns SSS-I and SSS-II can coexist. In (a-c) 
the combined state images (left), full power spectra (middle), 
and expanded power specta right) are presented at different 
temporal phases. The primary hexagonal symmetry (of mag- 
nitude fc<>) is broken by two different sets of wave vectors 
arranged in two 6-fold sublattices. The set corresponding to 
the SSS-I is enclosed in squares. The set corresponding to the 
SSS-II is enclosed in circles. This state was observed for 2/3 
driving and system parameters of v — 47 cS, cuo/(2n) — 25 
Hz and h = 0.31 cm, x = 62.5° and = 0°. 

This "cascaded" state appears as either a secondary 
bifurcation of an initial hexagon pattern or as a bifur- 
cation from a 12-fold quasipattern state. In its spatial 
power spectrum the two mechanisms that appeared in 
both SSS-I and SSS-II type patterns are cascaded and 
appear at different scales. An SSS-II type resonance oc- 
curs where wavevectors ql with magnitudes qi = K/y/3 
appear. A third vector Q, however, is also present. Q 
is half the magnitude of ql and, echoing the mechanism 
forming SSS-I states, breaks the 6-fold symmetry of the 
ql by aligning itself parallel to a single vector, q. 

Figure |l8| b demonstrates that this state has an overall 
temporal periodicity of 2t:/ujq. The two images (b left, 
right) were taken at an interval of tt/cuq- The images look 
exactly the same but with a transverse displacement of 
the d = 4tt/K length scale. The symmetry breaking by 
the Q is similar to the SSS-I mechanism where a mode 
with a phase of ir/2 produces a displacive effect, as shown 
in Figure |lqb. The peaks are of high amplitude and 
have the characteristic shape of the oscillons described 
by Arbell and Fineberg [|19| modified by the asymmetry 
that is produced by the displacive effect of the Q mode. 
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FIG. 18. (a) The image (left) and power spectrum (right) 
of "cascaded" SSS states appearing for 4/5 driving near Xc 
The two mechanisms appearing in SSS-I and SSS-II coexist 
at different scales. A hexagonal sublattice of qt — K/y3 bi- 
furcates from the original hexagonal pattern. In addition, a 
third wavevector Q ^ q breaks the symmetry of this sublat- 
tice in a single direction. This state has an overall temporal 
periodicity of 2tt/uio- This is demonstrated in (b) where two 
images (b left, right) were taken at an interval of ir/uio. The 
images are displaced by a d — 4ir/K length scale. The hori- 
zontal white line indicates the location of the peaks in the first 
phase (b left), (c) A the side view of this state. The peaks 
are of high amplitude and have the characteristic shape of the 
oscillons described by Arbell and Fineberg Jl9| with a small 
asymmetry in the direction of Q. These images were observed 
for system parameters of 75/60 Hz, v = 23 cS, h = 0.2 cm, 
X = 56° and = 0°. 

As demonstrated by Figure U9, different wavevectors 
are dominant at different temporal phases. These states 
were not observed in experiments with higher m/n values 
such as 6/7 and 8/9 driving ratio. 
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FIG. 19. Images (a-c top)and the corresponding power 
spectrum (a-c bottom) of different temporal phases of the SSS 
type pattern shown in Figure hq that appears in 4/5 driving 
experiment. The original hexagonal wave vector is drawn as 
the outer circle in the power spectrum shown (a,b bottom) . It 
can be seen that the different wave vectors have different rel- 
ative amplitudes at different temporal phases. In (a) the two 
smaller wave vectors Q and q have higher intensities than the 
wavevector K along the outer circle, where in (c) the center 
circle of radius q has a very small intensity. 



V. DOUBLE HEXAGON SUPERLATTICES (DHS) 

Historically, the first two-frequency experiment fo- 
cused on quasicrystalline patterns with 12-fold symme- 
try. These were observed in the vicinity of \c in systems 
driven with an even/odd driving ratio and were found 
to bifurcate either from the flat state or as a second 
bifurcation from the harmonic hexagon region. In this 
section we will describe both the 12-fold quasipatterns 
mentioned above and two other superlattice states (see 
Figure EOh . All of these states share a harmonic tempo- 
ral response for odd/even driving and appear in the same 
general area of phase space. An important factor in pat- 
tern selection is the nonlinear coupling coefficie nt tha t 
depends strongly on resonant locking (see Sec. ID2J ). 
In contrast to the SSS states, we will see below that no 
new wavenumbers are needed to construct these states. 
This class of superlattice patterns are differentiated by 
both the relative angular orientation and spatial phase 
relations of the critical wavevectors that form them. 

The 12-fold quasipattern (see Figure po|a) was ob- 
served for 4/5 driving in the vicinity of \c- This state's 
temporal behavior is harmonic with respect to luq This 
state can be formed near the linear threshold for small 
values of e ~ 0.01. Increasing the amplitude A causes a 
bifurcation to the cascaded-type superlattice described 
in IV C (see Figure O). The 12-fold quasipattern may 
also be understood as being formed by two sets of hexag- 
onal wavevector at 30° to each other. As in all quasicrys- 
tals, this state does not have long range order or a well 
defined sub-unit cell. It is interesting to note that unlike 
the quasipattern described in |17| that appeared only for 
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75° the quasipattern we describe exists for 



0°. 
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FIG. 20. When using even/odd driving, a class of tempo- 
rally harmonic states exists that is composed of two distinct 
hexagonal sublattices (delineated by squares and circles) of 
size K, the critical wavenumber of the harmonic driving com- 
ponent, (a) Double Hexagon State (DHS) formed by both the 
critical hexagonal lattice and an additional hexagonal sub- 
lattice of smaller amplitude which is rotated by an angle of 
~ 22°. (b) 12-fold quasicrystalline state (c) A similar ("SL1") 
state observed by Kudrolli et al. pl| for 6/7 driving. This 
state is formed by a resonance similar to (a) but with spa- 
tial phase of 2-7r/3 in each of the sublattice components, (a) 
and (b) were obtained for 4/5 (60/75 Hz) forcing, (c) was 
reproduced with permission from plj. 
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FIG. 21. Using our numerical modelling technique, we 
calculated the images of a perfect double hexagon state 
(DHS) formed by two hexagonal lattices rotated by 22°. (a) 
on — Pi = (see Eq. (H)) shows a negative crater state and 
(b) a positive peak state, (c) shows the calculated image of 
a 12-fold quasi-pattern. (d) A calculated DHS state with 
spatial phases of Oi = and /3i = 27r/3 as in Figure p0|(c). 

Another superlattice pattern (Figure ^0|c) which occurs 
in this regime of phase space was first studied by Kudrolli 
et al. with 6/7 driving ratio with loq/{2h) — lQAAHz, 
v = 2QCSt, and h — 0.3cm. They found a hexagonal su- 
perlattice pattern composed of two hexagonal sets of crit- 
ical wavevectors with harmonic temporal response. The 



equation that describes the surface height function of this 
pattern can be generally written down as: 



h(r) = Y^ A * coa{K t ■ r + a,) + ^ B. b cos(i?; ■ r + ft), 

(6) 



i=l 



1=1 



where \K\ = \K'\ = k c and 

and K[ can be obtained by rotating Ki by an angle of 
9 = 22°. In Figure pl|a,b,d we show simulated images of 
this equation for different values of j3 and Ai/Bi. 

The pattern (Figure |2C|c) described by Kudrolli et 
al. was found to have an angle 9 with the value 
9 = 2sin~ 1 (l/2v / 7) w 22°, equal amplitude coefficients 
1-^*1 = \Bi\ for i = 1, 2, 3 and spatial phase angles a, = 0° 
and ft = 120°. For this special value of 9, resonance con- 



ditions such as 2K 1 , 



K's 



2Ki — Kj, were shown |21 



to exist. 

Double hexagon states formed by the superposition 
of two hexagon sets of critical wavevectors oriented at 
a relative angle 9 = 2 sin" 1 (\j2\fl) are one example 
p3| , p0[ of a wider class of superlattices. This superlat- 
tice class is composed of a periodic lattice formed by a 
wavenumber smaller than the critical wavenumber of the 
excited surface waves. This smaller wavenumber corre- 
sponds to the six wavevectors formed by the difference 
between adjacent wavevectors of the two hexagonal sets 
i.e. by Qi=Ki-K[ for i = 1,2. ..6. Any two of these 
wavevectors, say Qi and Q 2 , are primitive vectors of the 
lattice and the vectors Ki and K[ are points on the Q~ 
lattice given by Ki — ri\Q\ + n 2 Q2- Silber and Proc- 
tor |23] show that only a discrete countable set of Qi 



where K/Q = yr 



n| — niri2 can satisfy this condi- 
tion. The angle 9 between the two lattices is given by 
9 = cos" 1 ( "/ 3 Kl " 2 , "jl ). In this formulation the ex- 
perimentally observed patterns are obtained for n\ = 3, 
n 2 = 2 and K/Q = \/7, giving 9 rs 22°. Only the sim- 
plest DHS have been observed to date. 

One characteristic aspect of the DHS states is the rel- 
ative spatial phase of the different wavevectors. The pat- 
tern described by Kudrolli et al. [Bl[ (coined "SL-1", see 
Figure |20p) consisted of triangular unit cells that were 
produced by a DHS with a t = 0° and ft = 120°. Silber 
and Proctor describe, in simulations of thermal convec- 
tion, a DHS state formed with both a* = 0° and ft = 0° 
fl23| . The DHS superlattice in our experiments with a 2/3 
driving ratio (see Figure |2(]a) differs from these super- 
lattices in that the two hexagonal sublattices that form 
them possess different amplitudes. States similar to these 
been observed in nonlinear optics |47fl . In addition, in our 
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case, there is no spatial phase difference between the two 
sets of hexagonal wavevectors (a* = (3j for all i,j). 

Although this state is stable close to the threshold, in- 
creasing the driving amplitude results in the appearance 
of many defects and eventual temporal disorder. At high 
driving amplitudes the spatial symmetry is hard to dis- 
cern due to the many defects and domains within the 
fluid cell. In Figure E3 we show a typical DHS state at 
different temporal phases. The state oscillates between 
a crater-like phase (a) and a positive oscillon-like phase 
with a frequency of uq. High amplitude oscillons appear 
in the center of the hexagonal sub-unit cell of the DHS 
state. Similar oscillons were also observed for the 12-fold 
quasipattern state shown in Figure EOp pj| . 



temporal resonance. 




22. Typical images of the double hexagon state (DHS) 
taken at different temporal phases, (a) shows a crater phase 
(b) an intermediate phase and (c) a peak phase for 2/3 driving 
with wo = 30 Hz, h = 1.55 cm, v = 23 cS and X = 60.4°. The 
time interval between (a) and (b) is 0.25T and between (a) 
to (c) is (0.4T) 

Let us summarize the common characteristics of DHS 
states. All of the patterns reveal a harmonic time depen- 
dence and the lack of any fundamental wavevectors other 
than the critical wavevectors excited by the harmonic fre- 
quency. All of these patterns were observed to be in the 
vicinity of \c on the harmonic side. In all cases, the mix- 
ing angle was greater than that needed for obtaining SSS 
states and within the range x° > x> (Xc~ 12°). Another 
common characteristic of these states is that all exist in 
the vicinity of a first-order transition of the hexagonal 
patterns from the featureless fluid state. This implies 
that quadratic interactions can play an important part 
in describing these states. Finally, two of these patterns 
also generate oscillons as described in detail in |19| ]. 



VI. TWO MODE SUPERLATTICES (2MS) 

In the vicinity of the critical mixing angle, Xci t w0 
modes with different wavelengths can be excited con- 
currently. These modes can interact in different ways 
to produce a variety of different patterns. In the next 
two sections we will describe two distinct types of two- 
mode states that are formed near \c- The first of these 
states are two-mode superlattices (2MS). These states 
are formed by the interaction of the two linearly excited 
modes with a third "slaved" mode that is selected via a 











FIG. 23. 2MS patterns consist of a nonlinear superposition 
of both eigenmodes, fci and k^, excited by the driving frequen- 
cies ui\ and L02- Although Wi and u>2 determine the lengths of 
fci and &2, their relative orientations are determined by the 
condition: hj, = &2 — fci where the wavenumber, k$, is deter- 
mined by the temporal resonance condition: us — u>2 — u)%. 
(a) 2MS patterns for 50/80 Hz driving where ki with square 
symmetry is dominant, (b) fo-dominant 2MS with hexag- 
onal symmetry for 45/60 Hz driving and (c) fci-dominant 
states having hexagonal symmetry for 40/60 Hz. Shown are 
the spatial spectra (center) with the resonant triads (right) 
&3 = &2 — k\ , highlighted. 

Figure p3| shows the three main types of 2MS states; ki 
dominated (with either square or hexagonal symmetries) 
and k\ dominated (with hexagonal symmetry). 2MS 
states exist in both harmonic and subharmonic regions 
of phase space in the vicinity of \c (see Figure [|) . In the 
region of phase space between the square 2MS and the 
hexagonal 2MS a spatially disordered "unlocked" state 
exists. Unlocked states, which are formed by the same 
wavenumbers that form the 2MS, have no well-defined 
spatial symmetries. As in the case of SSS states, the 
transition to 2MS from either square or flat states is non- 
hysteretic and occurs via propagating fronts. 

The 2MS are qualitatively different than SSS states. 
They result from spatial phase locking of both k\ and 
ki whereas the SSS states result from a resonance con- 
dition that is independent of the ki mode. 2MS states 
are the most general of the superlattice states described 
here. They are observed for all types of driving parities 
(odd/even, even/odd and odd/odd) and, as shown in Fig- 
ure E3, appear with either square or hexagonal symmetry. 
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A. 2MS Resonance condition 

2MS spatial spectra, as shown in Figure £3| are com- 
posed of peaks of length k\ and k 2 and their linear combi- 
nations. The strongest secondary peaks, indicated in the 
figure, are given by k% — k 2 — k\ where the magnitude of 
&3 is consistent with the linear value of k calculated for 
a single frequency excitation at the difference frequency 
UI3 = ui 2 — 0J1 ■ Our calculated value of k^ was obtained 
using the linear single frequency code of H] at threshold. 
The difference between the calculated and measured val- 
ues of £>3 varies between 5 — 20%. This shift between the 
measured and calculated values is constant for a given 
value of the difference frequency, W3 , and systematically 
decreases as W3 increases. We believe the shift to be the 
result either finite size effects in the cell or the fact that 
the ui(k) used is the linear dispersion relation for a fea- 
tureless state (not one with pre-existing waves). 

A subtle point in the interpretation of the power spec- 
tra of 2MS states, is the evaluation of the magnitude of 
the two wavevectors that appear. As can be seen in Fig- 
ure £4] the values of k\ and k 2 measured in the vicinity of 
Xc are significantly different than the values of k excited 
by single frequency excitation. The two-frequency linear 
stability code of Tuckerman [ [l3|]35| ], reproduces this ef- 
fect and agrees to within 1 — 2% with the experimentally 
measured values (see Figure |4|). This enables us to ac- 
curately calculate values of both k\ and k 2 in the vicinity 

of Xc- 

The resonant conditions stated above suggest that the 
orientation of the wavevectors building the 2MS is se- 
lected by nonlinear interactions that are resonant both 
in space and time. Thus, the temporal resonance con- 
dition dictates the spatial orientation of the vectors k\ 
and k 2 . Such 3-wave resonant interactions have been 
predicted to occur in non-linear interactions of surface 
waves p8| and are well known in the physics of plasmas. 
Similar states were observed as a result of nonlinear mix- 
ing of a multiple-mode optical beam pjl . The selection 
of &3 via the temporal resonance condition yielding U3 
is non-trivial and can not be accounted for by exper- 
imental artifacts such as possible nonlinearities in the 
imaging. This 3-wave resonance condition occurred for 
all frequency ratios tested. 

As in SSS, the dominant 2MS wavevector retains its 
initial symmetry, while the relative orientation of the 
other linearly excited wavevectors are determined by the 
above resonance condition. For odd/even driving hexag- 
onal (square) symmetry dominates for X < Xc {x > Xc)- 
Thus, square 2MS (Figure |23p.) bifurcate from the k 2 
square pattern that dominates the x > Xc region. Sim- 
ilarly, hexagonal 2MS (Figure p3|b) bifurcate from the 
k\ hexagonal pattern that dominates the x < Xc region. 
In 3/4 driving, hexagonal 2MS states are excited whose 
dominant scale is that of the larger wavenumber, k 2 (Fig- 



ure |23]c). It is known that in single-frequency driving 
experiments different symmetries can arise for different 
system parameters even when the temporal behavior is 
solely subharmonic. In two-frequency driving experi- 
ments the parity of the dominant frequency does not au- 
tomatically dictate the symmetry selected. For example, 
in 5/8 driving (50/80 Hz) 2MS states are observed with 
square symmetry in the temporally harmonic region. 
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FIG. 24. A comparison of k (in cm 1 ), calculated using the 
two-frequency linear calculation [|13| (solid line) ,with mea- 
sured values of k as a function of \- Note that both values 
of k are nearly constant away from y c = 58°. Near \c, ki 
(corresponding to u> 1 = 40 Hz) decreases by nearly 10% while 
&2 (corresponding to u>2 = 60 Hz) increases slightly. 



B. Temporal behavior of 2MS 

Using our strobed imaging technique, we studied the 
spatio-temporal behavior of the patterns by inspecting 
them at different temporal phases. The 2MS inherit 
the basic temporal periodicity of both primary eigen- 
functions, which are periodic in time with a basic fre- 
quency of cither luq (harmonic) or luq/2 (subharmonic). 
The superposition of the two modes always displays a 
subharmonic response (ojq/2). Like the SSS, the state's 
appearance changes qualitatively with time, with the 
spatial wavenumbers having fixed locations. 2MS time- 
dependence stems from the temporal dependence of the 
amplitudes of these modes. In Figure Eq and Figure 26 



we 

see typical sequences of both square and hexagonal 2MS 
states taken at constant values of the driving parame- 
ters. The relative intensities of the different wavevec- 
tors change in the different phases; in Figure |25| e the k\ 
and &3 wavevectors are almost absent, resulting in a pat- 
tern whose dominant square symmetry has a 2ir/k 2 scale, 
while in the remainder of the phases all vectors appear. 
In Figure E0 the resonance condition allows k\ and k 2 
to be nearly co-linear and two sets of 6-fold wavevectors 
appear in the power spectrum. As in the square 2MS, one 
can see phases with significantly stronger k\ peaks (Fig- 
ure |26|(a)) as well as those where either the kj, (Figure 
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|26|(b)) or fc 2 (Figure ^6|(c)) are stronger. 

The relative stability of 2MS hexagonal states (e.g. 
Figure p3|b,c) is dependent on whether ki or fc 2 is dom- 
inant. Let us first consider hexagonal 2MS states where 
k\ dominates. This state is found for even/odd driving 
for x < Xc (see FigureH) and a typical time sequence 
is presented in Figure |27| . Again, different wavenumbers 
are dominant in different temporal phases. 








FIG. 25. A typical temporal sequence of square 2MS states 
taken at constant parameter values for 50/70 Hz driving, fci 
(middle circle), fe (outer circle) and fe (inner circle) are in- 
dicated in (a). In (d) the k$ vectors are absent but enclosed 
in squares are vectors that are the vector sum k' 3 + ft3. Sys- 
tem parameters are ujo/(2tt) = 20 Hz, v = 23 cS,</> = 0°, 
X = 65.4°, and h = 0.155 cm. 




FIG. 26. A typical sequence of hexagonal 2MS (a-e) at dif- 
ferent temporal phases taken at constant parameter values 
for 60/80 Hz driving. Circles of radii fci (middle) , &2 (outer) 
and &3 (inner) are drawn in (c). In (b) the &3 and fci vectors 
interact to produce the vectors enclosed squares. k\ and hi 
form two sets of hexagonally arranged vectors. 

The region of phase space where this state appears is not 
as large as the fc 2 dominant 2MS, but these states ap- 
pear for all even/odd driving ratios used. As in the case 
of square 2MS and the 2kR states discussed in Section 



VII, these states are not observed for all combinations 
of 4> and S/h used. It is difficult to obtain a hexago- 
nal 2MS state of this kind that extends over the entire 
system and fci-dominant hexagonal 2MS states generally 
occur within domains. Thus the spatial spectra (as seen 
in Figure 27) appears sometimes smeared, as the various 
domains have different angular orientations. Although 



the scenario described in Figure £7] is typical, we have 
observed stable global k\ -dominated 2MS state for par- 
ticular values of 4> and 8/h (see e.g. Figure p2| ). 

A typical time sequence and power spectra for the ki 
dominant hexagonal 2MS state (shown in Figure E3c) is 
presented in Figure Era. Here, in contrast to k\ dominant 
hexagonal 2MS, the pattern is global with a well-defined 
symmetry in all of its temporal phases. Hexagonal sym- 
metry was not seen to be preferred for all odd/even driv- 
ing and was only observed for the simple ratios of 1/2 and 
3/4. This might be a feature of simple driving ratios that 
can sometimes have unique properties which are related 
to temporal locking ( [[ll|]33; l). 
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FIG. 27. A typical sequence of fci-dominant hexagonal 
2MS states taken at constant parameter values for 40/70 
Hz driving. These states generally appear in domains. 
(Wavenumbers are noted by arrows). 

2MS states are not always as highly ordered as those 
in Figure |25| and Figure g6|. When highly ordered 2MS 
states occur, the resonance condition for these parame- 
ters allows ki ki and k% to be nearly co-linear. When the 
angles between the resonant wavevectors are not small 
the pattern is usually not global and a number of do- 
mains coexist. In highly ordered states we find a locking 
of the wavevectors' magnitude to small natural number 
ratios. The locking ratios of k\ : fc 2 : k$ — 1 : 3 : 4 and 
1:4:5 were, respectively, obtained in Figure pojand Figure 
EO. This "spatial locking" seems to stabilize global pat- 
terns in a way similar to the global SSS states presented 



in Section III. 



The local laser probe technique provides more detailed, 
quantitative, information of the temporal behavior of 
2MS states. In Figure g8| we show one component of 
the slope of surface waves obtained for 4/5 driving. In 
Figure |28|(a) we show the typical waveform of a hexag- 
onal pattern for x < Xc- Besides the strong component 
at Wi/2 = 4^o/2 and its harmonics (40,60,80 Hz ...), 
peaks appear at values of jujo,j — 1,2.... Those peaks 
are expected from the linear theory. In Figure |28|(b) we 
see a typical waveform of a square pattern for x > Xc- 
Although the strongest frequency response is at cl>2/2, 
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the basic frequency, ujq/2, together with its higher har- 
monics also appear. Linear theory predicts a different 
distribution of energy in the peaks for the unstable mode 
at threshold. Zhang and Vinals |1T|] non-linear theory 
accounts for the peak's strength in a semi-quantitative 
way, as it correctly predicts the ordering of the strongest 
peaks. 
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FIG. 28. The time dependence of a two-frequency exper- 
iment with even/odd driving was studied using the reflec- 
tion of a laser by the surface waves. The temporal response 
is subharmonic (harmonic) with respect to ujq = 10 Hz for 
X < Xc (x > Xc)- Time series (left) and power spectra 
(right) for three different regimes: (a) Typical temporal re- 
sponse for hexagonal pattern found for X < Xc - Although 
the 4a>o/2 = 20 Hz component is the strongest, an ojo = 10 
Hz harmonic component with all of its higher harmonics is 
present, (b) Typical time dependence for square patterns 
(found for X > Xc)- The strongest response is for 4wo/2 = 25 
Hz as expected. A weaker response at the basic subharmonic 
Wo/2 = 5 Hz is observed together with stronger responses 
at other harmonics. The ordering of the harmonics agrees 
qualitatively with jllj. (c) Typical temporal response of a 
square 2MS state at Xc . The power spectrum shows a strong 
response at both uj\ (20 Hz) and u>2 (25 Hz). All other har- 
monics are present with a different power distribution than in 
(a) and (b). System parameters here are: 50/40-ffz, v = 23 
cS, h = 0.155 cm, (f> = 0° and X = 40° (a), X = 60° (b) and 
Xc = 53.5° (c). 



In Figure |28|(c) the temporal response of the 2MS state 
is presented. This state has a temporal response that 
includes both the frequencies of the harmonic mode (a) 
for x < Xc and those of the subharmonic one (b) for 
X > Xc- Although the power spectrum shows that both 
main peaks at uq/2 and tu 2 /2 are of the same strength, 
their relative strengths can vary with the location of the 
laser probe. 

We have seen that the 2MS state contains the two 



linearly excited eigenmodes both in space and in time. 
Does each mode keep its distinct space-time behavior, 
or is there a complete mixing of the spatial and tempo- 
ral components? To clarify this question, we write down 
a simple model for the surface height of a square 2MS 
state. For simplicity, we will assume a square 2MS state 
where both k\ and k 2 are collinear with the same spatial 
phase: 

(7) 



h(r,t) =.Fi(i)(cos(feia:) 

+F 2 (t)(cos(k 2 x) 



-cos(feiy)) 

- cos(k 2 yj) 



dh 

— = fei-Fi (t) cos(kix) - 
ox 

— = kiFi(t) cos(kiy) 
ay 



k 2 F 2 (t) cos(k 2 x) 
k 2 F 2 (t) cos(k 2 y) 



where iq(t) ~ cos(wii/2) + ... and F 2 (t) ~ (cos(u 2 t/2) + 
.... It is easily seen that if one selects a point (x.y) = 
(ir/2k 2 ,Tr/2ki) then each of the components of the par- 
tial derivative dh/dx, dh/dy has a temporal dependence 
of iq(i) and F 2 (t) respectively. Both Fi(t) and F 2 (t) in- 
clude time dependent terms retaining the parity of the 
linear modes: odd multiples of (p+ l/2)wo for a subhar- 
monic response and even multiples of pu>o for harmonic 
response. 
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FIG. 29. The time dependence of a square 2MS state mea- 
sured via the reflection of a laser at the fluid surface. Time 
series (left) of surface gradients, d y h (a) and d x h (b) are shown 
with their power spectra (right), (c) d y h as a function of d x h. 
The separation of the temporal response (a,b) was obtained 
by choosing x and y along the symmetry directions. d x h is 
dominated by a 3tJo/2 = 15 Hz response whereas d y h under- 
goes a 5o;q/2 = 25 Hz response. System parameters used were 
50/30 Hz, v = 23 cS, h = 0.155 cm, <f) = 0° and X c = 56.3°. 

As suggested by Eq. R, experiments (see Figure |29| ) 
confirm the separation of the time dependence of the 
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two modes within the 2MS state. The d x h component 
has a strong u>i/2 response whereas the d y h component's 
strongest peak is at W2/2. This strong separation of spa- 
tial time dependencies can only be observed for a few 
points in the (x,y) plan and, in general, the two frequen- 
cies are mixed. 



C. The "unlocked" state and transition regions 

Let us now consider the "unlocked" state that appears 
in the near vicinity of Xc- In Figure |3fj we present a 
typical time series of the "unlocked" state and its cor- 
responding spatial spectra. The spatial behavior of the 
state varies rapidly, over time scales of order (2tt/u>o). I n 
contrast to the SSS and 2MS states, in the "unlocked" 
state no orientational order is apparent. 
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FIG. 30. In the "unlocked" state no orientational order is 
apparent. As in 2MS states, k\ , &2 and kz exist simulta- 
neously but "unlocked" state spectra are diffuse and show 
little angular correlation, (top) Typical views of an unlocked 
state with corresponding spatial spectra (bottom) observed 
for 40/60 Hz driving at different times and fixed driving pa- 
rameters. The different values of k appear with different mag- 
nitudes and orientations in each phase. Their lengths are in- 
dicated by the arrows. 

Both fci and fe exist simultaneously in their spatial spec- 
tra but spatial mode locking does not occur as in the 2MS 
states. This is evident in their power spectra, where, gen- 
erally, entire circles of radii k\ and fe appear. As Figure 
pfj| d indicates, additional peaks of wavenumber k% corre- 
sponding to L03 = u>2 — u>\ are sometimes observed. The 
"unlocked state" is a well-defined state that exists in a 
relatively wide region of phase space. This can be seen 
by defining (as in Blnl) the following "orientational cor- 
relation function", Cf.(6) for each value of k: 

C k (0) = E a [/ fc (o) ■ f k (a + 0)] /£«[/*(«) • /*(«)] (8) 

where fk (a) is the Fourier transform of the wavenumber 
k at the polar angle a. The correlation amplitude, Qk = 
l/2[maxCk(9) — minCk(6)], varies between and 1 for, 
respectively, minimal and maximal orientational order. 



As Figure |3l] shows, both the 2MS and hexagonal states 
have clear orientational order while very little residual 
order is apparent in the unlocked state. 

As is apparent from Figure pU the orientational am- 
plitude drops sharply as the boundary between the 2MS 
and "unlocked" phases is crossed. This is demonstrated 
in Figure p2|a-d which corresponds to the range of 5.1- 
5.25g in Figure Bl]. Between the pure hexagonal state 
and the "unlocked" state the hexagonal 2MS state ex- 
ists. This state displays strong k\ dominant hexagonal 
symmetry at most temporal phases but at other temporal 
phases some of the power is in the ki mode. Since Qk is 
an average of equally time separated temporal phases the 
hexagonal 2MS state has a Qk which is slightly smaller 
than that of the pure hexagon state. 



0.75- 



| 
a, 



'f 



0.50 



■•r 



Q O.25 



.: ! "Vv' ! 



4.5 5.0- 5.5 7 

/ Driving Amplitude; (g) 



6.0 



-.-$&*£•. 

-:•*£!$£;* 




90 ISO 90 180 
UNLOCKED HEXAGON 



FIG. 31. Within the "unlocked" states, the correlation am- 
plitude, Qk drops significantly and the angular correlation 
function, Ck{0), indicates only a small amount of residual or- 
der, (upper) Qk 1 (triangles) and Qk 2 (squares), averaged over 
a single period along the line \ — 58° for 40/60 Hz driving as 
a function of the driving amplitude, A. This line traverses the 
square, square-2MS, "unlocked", and hexagonal-2MS phases. 
Typical patterns in each phase are displayed (center). The 
symmetry of the different phases is highlighted by Ck (0) for 
fci (grey line) and &2 (black line) computed for typical states 
{0 is in degrees) . The power (log scale) of k-z relative to k\ in 
each C h (0) is 8.3 (square), 1.1 (2MS), 0.9 ("unlocked") and 
0.2 (hexagon). 

The parameter, Qk, does not completely characterize 
the different transitions, since spatial FFT power spec- 
tra cannot differentiate between global and local order- 
ing. A closer look at real-space images of the transitions 
presented in Figure 31 that correspond to square-2MS, 
"unlocked" and Hexagonal-2MS are shown respectively 
in Figure B2^ and Figure p3. Both transitions involve an 
advancing front that separates two well-defined domains. 
The clear separation between unlocked and 2MS domains 
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provides further evidence that the unlocked state is in- 
deed a distinct nonlinear state and not, simply, a transi- 
tion region. 

The transitions between unlocked states and the 2MS 
states with different symmetries differ in two ways. The 
sensitivity to any change of the driving parameters is 
much higher for the unlocked to 2MS-Hexagon transi- 
tion. Whereas the transition between square 2MS and 
"unlocked" states occurs for a relative change of ampli- 
tude of < 5%, the transition between hexagonal 2MS and 
hexagons can occur via a change smaller than 0.25%. The 
time scales of the induced transitions are also different. 
The first transition takes place in a nearly quasi-static 
reversible way, whereas the second transition (as shown 
in Figure p3| ) can occur over typical time scales of 50- 
1000 oscillation periods with a hysteresis of less than 
0.1%. This rather sharp transition is, perhaps, due to 
the effects of the quadratic interactions inherent in the 
harmonic states. The precise duration of this transition 
depends on the initial and final driving parameters. As 
is typical of front propagation processes, the deeper one 
is within the hexagonal regime the faster the transition 
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FIG. 32. The transition between the square 2MS state 
and the "unlocked" state occurs through a gradual process 
in which both states co-exist in different domains. At the 
transition, increasing the amplitude for constant \ constant 
results in (a) the global 2MS state is first disturbed by small 
defects at the cell's rim. (b,c) The disturbance spreads to the 
cell's center until, finally, the entire pattern is in "unlocked" 
state(d). This process can also occur in the reverse direction. 



resonant mechanisms found in the vicinity of \<P- 

In this section we will describe the spontaneous forma- 
tion of nonlinear rhomboidal patterns, formed as a re- 
sult of solely spatially resonant 3-wave coupling between 
wavevectors with distinctly different wavenumbers. This 
state has only previously been observed in a nonlinear 
optical system where the orientations of the interacting 



wavevectors 51 were externally imposed. Rhomboidal 
patterns have also been recently observed in paramet- 
rically driven ferrofluids |45|. The rhomboidal patterns 
described below spontaneously couple the two circles of 
linearly degenerate wavenumbers. These states qualita- 
tively differ from 2MS states in that they are composed 
solely of the linearly excited wavevectors, k\ and &2, in 
contrast to the additional slaved mode (k^) necessary for 
2MS formation. 

The rhomboidal states observed in this system re- 
sult from the non-linear interaction of k\ and &2, 
which are waves with significantly different wavelengths. 
Such states have been observed numerically in a Swift - 
Hohenberg type models |25|,|27J as discussed in section Q. 
They have also been anticipated in anisotropic models 
where two degenerate wavevectors are resonant with an 
externally imposed wavenumber |B2| , in nonlinear optical 
systems J5^j54[] , and in the analysis of the Faraday insta- 
bility excited with two frequencies p3|,B0||. Both rhom- 
boidal states and superlattice patterns have also been 
recently predicted to occur as a result of two bistable 
modes coupling to a zero-mode |55| . 

The rhomboidal states observed in this system differ 
distinctly from rhomboids resulting from slightly "dis- 
torted" hexagonal states J5(3|)7|]. "Distorted hexagons", 
predicted to be stable in models with derivative-coupled 
quadratic terms, may arise due to either initial or bound- 
ary conditions |58| . These states have been observed in 
reaction-diffusion systems |57| ] , convection in an imposed 
shear flow |59|, and flux line lattices in superconductors 




iij: : : : :. : :. : !. : i:3 

(a)'""-" W(hT* F (c)^ ^(d)^ 

FIG. 33. A time sequence of the transition from the reg- 
ular hexagonal state (a) to Hexagonal 2MS states is shown. 
The hi (60 Hz) component propagates via a front (b,c ) un- 
til it dominates the entire plate and co— exists with the large 
length scale (d) . The driving parameters are constant for this 
sequence. 



VII. RHOMBOIDAL STATES (2KR) AND 
QUASI-PATTERNS 

In the previous section we studied the three wave res- 
onance mechanism that can form 2MS states from the 
interaction: ^3 = k% — k\ and W3 — 102 — ^1 • Is the only 



A. Experimental conditions 

We have observed rhomboid states using both Dow- 
Corning 200 oils with kinematic viscosities, v, of 8.7, 23, 
47, 87 cS (at 30° C) and TKO-77 vacuum pump fluid with 
v = 184 cS (at 33°C) in fluid layers whose depth varied 
between 0.1 < h < 0.55 cm. The 2 wavenumber rhombic 
(2kR) states described in this section were all generated 
with m/n = 2/3 and 12 < ojq/2tt < 45 Hz or m/n = 4/5 
with 10 < ujq/2tt < 20 Hz. Frequency combinations of 
25/50, 30/60, 40/80, 40/56, 45/63, 48/68, 50/70, 55/75, 
55/77, 30/50, 40/70, 45/60, 50/65, 50/80, 52/68, 60/84, 
60/100 did not generate 2kR states. 
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transition to DHS and oscillon states also occurs as in 
Figure ||[ 
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FIG. 34. (a) A typical 2k Rhomboid (2kR) state observed 
for 2/3 driving and u>o/(2tt) = 25 Hz, v = 23 cS, and h — 0.2 
cm. The axis y is defined by the direction of ki. (b) The 
power spectrum (c) illustrates the simple resonance condition 
k' 2 — k2 = fci that generates these states. k' 2 and k 2 correspond 
to u>2 = Suio/2 while ki corresponds to <jj\ — 2uoo/2. 9 here is 
41°. 

In Figure p3 we present a typical image of a 2kR 
pattern. Although k\ and k 2 are the linearly unstable 
wavenumbers E3] excited by u>\ and cu 2 , their orienta- 
tion is determined by the 3-wave non-linear interaction 
that yields the resonant triad: k' 2 — k 2 — fci, where 
\k' 2 \ = \k 2 \ = k 2 . The higher harmonics in the figure 
may either be real or could occur as an artifact of the 
imaging. An additive three-wave resonance occurs for 
1/2 forcing where the resonance k[ + k\ — k 2 governs 
the selected patt ern. In this case, the resulting pattern 
(see section VIII ) is a superposition of hexagonal lattices 
composed of the two scales. 

The observation of 2kR states solely for driving ratios 
2/3, 4/5 and 1/2 is entirely c onsist ent with Silber and 
Skeldon's |3(j] predictions (see [D2| ) that 3-wave inter- 
actions coupling the wavenumbers k\ and k 2 are only al- 
lowed when two odd-parity waves are coupled to a wave 
with even parity. Thus, k! x + k\ = k 2 coupling is al- 
lowed for odd/even frequency ratios such as 1/2, and 
k 2 — k 2 = k\ coupling occurs in even/odd forcing, such 
as 2/3 and 4/5. It is interesting that we have not ob- 
served these states for other frequency ratios. It is possi- 
ble that 2MS states are preferred for all but the simplest 
frequency ratios since, for higher ratios, linearly stable 
tongues corresponding to wavenumbers close to the value 
of &3 (defined by ui^ — oj 2 — wi) are more dense. 

A typical phase space in which 2kR states are observed 
is presented in Figure [| (right). For values of x that are 
far from Xc the phase diagram is similar those described 
in sections III and VI. The 2kR state exists in the near 



vicinity of Xc and replaces both the 2MS and unlocked 
states. This region is bounded for x > Xc by squares 
and for x < Xc by k\ -dominant DHS states that are 
mixed with oscillons pffl (see [[L9[). At higher values of 
v and h (e.g. v = 47 cS, h = 0.3 cm, v = 87 cS, h = 
0.5 cm) square and hexagonal patterns only exist near 
onset. Upon increase of A both types of patterns become 
rolls. The 2kR state is, however, unaffected by the state 
preceding it. They appear for a similar range of <p. The 
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FIG. 35. A common phenomenon observed in the 2kR re- 
gion of phase space is the formation of two domains of 2kR 
states (left) which are oriented with a relative angle of 9 that 
is also the angle between the two larger vectors k 2 , k 2 forming 
the 2kR state. A similar phenomenon can also be observed in 
domains formed by triblock copolymers |6l|. 

2kR states are not always correlated throughout the 
entire system. At relatively low values of A, 2kR states 
can sometimes be found in 2 or 3 domains, as shown in 
Figure pq. The angle separating two such domains is 
identical to the angle 9 between k 2 and k' 2 , as defined 
in Figure p4|(b) . This type of domain separation is also 
observed in "knitting patterns" pi]], formed by triblock 
copolymers near a bistable point. (These new materials 
have a reciprocal lattice structure similar to 2kR states). 




FIG. 36. The transition from the 2kR state to the DHS 
state occurs through an "explosive" process in which a spatial 
domain of hexagons forms and rapidly increases in size. Six 
images (a-f) of this process taken at 1.5 second intervals for 
an experiment performed with 66/44 Hz driving for v — 23 
cS, and h = 0.2 cm at mixing angle \ — 70°. Note areas of 
mixed 2kR and hexagonal domains around the perimeter of 
the hexagonal area (clearly seen in upper half of (e)). In this 
process the driving parameters are fixed. 

Let us return to the description of the typical phase 
space. As increasing the driving amplitude A strengthens 
the nonlinear interactions between the waves, 2kR do- 
mains coalesce at higher levels of excitation. Further in- 
crease of A yields a hysteretic bifurcation to the "Double 
Hexagonal Superlattice" state (DHS) , where oscillons are 
formed at the maxima of the pattern (see jl9[ ) . In Figure 
p6| we present a time sequence of the transition from a 
global 2kR state to a DHS state by means of rapid front 
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propagation. Only at the final stages of this process does 
the pattern bifurcate into the high-amplitude state con- 
sisting of oscillons superimposed on a DHS background 
(see Sec. |y| and ]19|]). The opposite transition from the 
DHS to the 2kR state has a qualitatively different charac- 
ter. Small A^-dominated wavelength domains penetrate 
the DHS from the perimeter of the cell in a way similar 
to the transition from hexagons to SSS. The transition 
is not reversible and can have a small hysteresis (under 
1%). 

Both 2MS and 2kR states are observed for the driving 
ratios 2/3 and 4/5. In p2f the dimensionless parame- 
ter 5/h = (v/uave) 1 / 2 /h (where ui ave = (uj 2 - Wi)/2) 
was shown to govern the selection between the two pat- 
terns. For = 0, 2MS/ "Unlocked" states exist above 
5/h ~ 0.12 - 0.17, while 2kR states exist below. This 
critical range of 5/h was obtained for a broad range of 
both h (0.1 < h < 0.6cm) and v (8.7 < v < 186cS). 
The parameter 5/h is the ratio of two important physi- 
cal scales of the system: the ratio of the viscous bound- 
ary layer length, where the flow is rotational, to the fluid 
layer's height. This parameter, in essence, defines the 
region of applicability of the Zhang and Vinals, quasipo- 
tential approximation. (5/h) 2 is the ratio of the dissi- 
pative time to the driving time scales. In Lioubashevski 
et al. this number was critical for determining pattern 
selection by single frequency excitations |j^,[56). Thus, 
this transition suggests that high dissipation in the sys- 
tem favors the 2MS over the 2kR states. This may result 
from the (linear) broadening of unstable tongues |Q that 
occurs when dissipation in the system is increased. This 
broadening would make the linearly stable wavenumber 
observed in 2MS states, k 3 , more accessible. 
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FIG. 37. When 5/h is at the transition between 2kR and 
2MS stability zones (b) 2kR and (c) 2MS exist in the same 
vicinity. (A = 6.59g, \ = 60.6° for 2kR and A = 6.53g, 
X — 61.4° for 2MS). (a) A transient state combining the two 
resonances was obtained for the parameters of (b). Parame- 
ters here were 2:3 driving, luo — 25 Hz, and 5/h — 0.12. 

In the transition regime of 5/h, both 2MS and 2kR 
states can coexist, as shown in Figure B7Uor 5/h = 0.12. 



In this region, at fixed values of 5/h, small changes in 
either A or x can result in globally stable states of either 
type. As shown in Figure |37Ja, transient states in which 
both states are present can also result in this regime. 
In these states both resonant mechanisms can operate 
concurrently in different spatial regions of the fluid cell. 
As the two driving frequencies are commensurate, the 
phase variable, <j> in Eq. YU is a relevant control parameter. 
As was shown in the experimental work of Muller and 
Edwards and Fauve (see Sec. FG) changing <f> can affect 
pattern selection. Typically, the 2kR state exists over 
the range —20° < cj> < +15°. The phase space presented 
in Figure [| (right) is typical for 0.16 < h < 0.22 cm and 
v = 23 cS. In Figure |38| we show how changing the angle 
(j> causes the pattern to change from the 2kR to the 2MS 
state. The transition is not abrupt and in some regions 
localized patches of both states can co-exist. In general, 
the size of the region in phase space where a single 2k R 
domain exists decreases with the distance from (j> = 0° . 
Because <j> = 180° is equivalent, for 2/3 driving, to <j> = 
0°, the 2kR is also stable at this angle. The strong effect 
of changes in <p is consistent with the predictions of Zhang 
et al. (see Sec. ID1) who showed how <f> can affect 
the mode coupling function (3(9). Silber et al. |33|] have 
recently demonstrated that changing <f> affects the 2kR 
resonance by varying the non-linear coefficients of the 
model equations describing the system. 




FIG. 38. The effect of changing (f> in the driving can be seen 
in the sequence of images (left) and their power spectra (right) 
for different values of <f). While for = 0° (a) a 2kR state is 
stable and global, increasing cf> from (b) 20° and (c) 30° trans- 
forms the system to a stable 2MS state at (d) 90° . Increasing 
the angle to (e) 120°, (f) 150° reverses the process and fi- 
nally for <f) — 180° a global 2kR re-emerges. The 3 concentric 
circles indicate the magnitude of the three 2MS wavevectors 
(&i,fa, kz), while the gray arrows indicate the wave vectors 
of the 2kR state k\,k2- These measurements were performed 
for 2/3 driving with system parameters: ujo/(2ti) — 25 Hz, 
v = 23 cS, and h = 0.2 cm at mixing angle x = 70.5° 



26 



B. Temporal behavior 

The temporal behavior of the 2kR state is similar to 
that of the even/odd 2MS states, where the time de- 
pendence in different directions is qualitatively different. 
Like the even/odd 2MS states, different directions can 
have dominant subharmonic or harmonic components. In 



Figure 39 we see how the harmonic nature of the 2kR is 
manifested. A time translation of tt/uq shifts the ob- 
served pattern by a spatial translation of Tr/\k 2 + k' 2 \ in 
the x direction (as defined in Figure p4|a). The overall 
spatio-temporal behavior of the 2kR state is consistent 
with Eq. (§). 

h (t, x, y) — [a\ cos(cjoO + »2 cos{2ujQt) + ...] cos(/ci • x) 

+ [61 cos(oj t/2) + b 2 cos(3cj V 2 ) + •■■] x (9) 
[cos(fc 2 • y) + cos(k 2 ■ y)) . 

The fact that the two directions, x\\ki and y _L x (as 
defined in Figure p4T) each exhibit different time depen- 
dence is demonstrated in Figure |4C|. This typical time 
series of the x and y components of the surface gradient 
of this state at a single point is similar to that presented 
in Figure £9] for the 2MS state. By our choice of axes, 
the d x h component contains mainly the to\/2, 2uji/2 . . . 
peaks, while the dominant frequencies in the y direction 
are (lo 2 — lui)/2, lo 2 /2 . . .. Arbitrary x and y directions 
will contain both k± or k 2 eigenmodes. This orientational 
dependence may prove to be a general characteristic of 
superlattice states. When a single mode is dominant, the 
predictions of Zhang and Vinals pi[ are in good quanti- 
tative agreement with our measurements of the relative 
peak intensities. In the vicinity of Xc, when 2 modes are 
concurrently excited, this analysis does not apply and a 
new theoretical framework is needed. 
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FIG. 40. The time dependence of the 2kR state was studied 

using the reflection of a laser by a point on the fluid surface. 
The directional derivatives d x h and d y h as a function of time 
(left) and their corresponding power spectra (right) show that 
each direction has a different temporal dependence. (The x, y 
directions are defined in Figure B4a). d v h is dominated by the 
hi and k' 2 components with dominant frequency W2/2 and d x h 
corresponds to the k\ component with dominant frequency 
wi/2. This figure was taken from |22|]. The parameters for 
the above were u>q/(2tt) = 22 Hz, v = 23 cS, cj> = 0° and 
h = 0.2 cm 

In low viscosity fluids (y = 8.7 cS, 0.1 < h < 0.2 
cm, and S/h < 0.13) an interesting variant of the pure 
2kR state is observed whose symmetry changes with its 
temporal phase (see pSq]). At different temporal phases 
pure hexagonal, mixed hexagonal and 2kR phases can 
be seen. Interestingly, this state exists for a significantly 
broader range of <f> (—70° < <j> < 70°) than the pure 2kR 
state. Together with the vector triad characteristic of 
2kR states, coupling with the difference vector, ki — fei, 
is also observed in the spatial spectra of these states. 



C. Tuning of the resonant angles and quasi— pattern 
formation 





FIG. 39. A typical 2kR state taken for u /(2-k) = 25 Hz, 
v = 23 cS and h — 0.2 cm, in temporal phases separated by 
half of the common excitation time T/2 = 7r/wo- The time 
translation is equivalent to a spatial translation of magnitude 
d — 7r/(|/s2 + k' 2 \). This feature can be reproduced by Eq. H. 



The angle 6 between the two wavevectors k 2 and fc 2 , 
can be tuned by changing the different system parame- 
ters. Since the values of k\ and k 2 are roughly determined 
by the dispersion relation u)(k), the angle 9 can be varied 
by changing ojq, v, or h while leaving m/n constant. This 
is demonstrated in Figure |4l]. The value of 8 can be cal- 
culated for a specific mixing angle \ using the numerical 
method developed by Kumar and Tuckerman fllpf. 
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FIG. 41. The resonant angle between fc 2 an d &2 varies with 
the parameters of the experiment. We obtain angles of (a) 
36° for system parameters of ujo/{2n) = 25 Hz, v = 47 
cS, and h, = 0.25 cm, (b) 33° for system parameters of 
ujq/{2tt) — 15Hz, v = 87 cS, and h = 0.5 cm, (c) 32° for sys- 
tem parameters of wo/(27r) = 25 Hz, v — 87 cS, and h — 0.5 
cm and (d) 29° for system parameters of u>q/(2k) = 20 Hz, 
v = 184 cS, and h = 0.54 cm. All images were taken from a 
9x9 cm 2 square in the center of the circular plate. 





FIG. 42. (a) A 10-fold quasiperiodic pattern and its power 
spectrum is observed for uio/(2n) = 30 Hz, v — 47 cS and 
h — 0.33 cm. Circles are drawn with radii fei (inner) and 
&2 (outer). For these parameters 9 — 360°/ 10 and five 2kR 
regions combine to form the quasipattern. For ujo/{2it) = 30 
Hz, v — 23 cS and h = 0.2 cm (b) an nearly 8-fold symmetric 
pattern is seen. This pattern, however, is actually a deformed 
quasipattern since the 2kR value here of 9 — 41° does not 
evenly divide 360°. Instead of a single angle, two different 
angles 9\ = 41° and #2 = 49° are observed in the power 
spectrum. 

As suggested by J2J^(|, when rhomboidal states ex- 
ist the tuning of 9 can provide a well-defined mechanism 
in which quasipatterns of any desired symmetry may be 
formed. When 360/0 = p, an integer number of adja- 
cent triads can be formed. As conjugate pairs of triads 
are always formed, the integer p must always be even 
(as observed in |l^,gj). This is demonstrated in Fig- 
ure [42] where the formation of perfect 10-fold quasicrys- 
talline patterns and approximate 8-fold quasi-patterns 
occurs for values of 9 that are tuned to these resonant 
(360/0 = n) angles. As the power spectra of these states 
indicate, each of the inner circle of peaks of magnitude 
k\ is coupled by a triad resonance with two peaks of 
magnitude k 2 along the outer circle. When 9 rs 41°, a 
symmetric quasi-pattern is not possible and a distorted 
8-fold quasi-pattern, as shown in Figure Kt2b, occurs. 



D. Three frequency driving 

Having observed the distorted 8-fold quasipatterns de- 
scribed in Figure |42|b, we attempted to stabilize these 
"asymmetric" quasi-patterns by modifying the driving. 
Mullcr in |12J added a third frequency perturbation to 
break the spatial phase symmetry in the subharmonic 
regime and thereby control the transition between tri- 
angles and hexagons. This motivated us to add a third 
frequency in order to enable the excited wavevectors k\ 
and k 2 to spatially lock to the value of 9 = 45° for which 
8-fold quasi-patterns can naturally form. We used the 
following driving function: 

A[aicos(piui a t) + a 2 cos(p 2 uj t + <j>\) + a 3 cos(p 3 uj a t + (f> 2 )] 

(10) 

where the total driving amplitude is given by A and the 
normalized amplitude ratios by a\ : a 2 : 03 with a\ + 
a>2 + 0,3 = 1. Pi '■ P2 '■ P3 are the three frequency ratios 
Pi < Pi < P3 and 0i, fa are the phase differences with 
respect to the p\ components. 





FIG. 43. Two temporal phases of an 8-fold quasi-pattern 
observed for three-frequency driving (a,b). This state is tem- 
porally subharmonic and was observed in the region where the 
distorted 8-fold quasi-patterns were observed (see Figure kt3. 
This state was observed in a 50/75/100 Hz experiment with 
driving amplitude ratio a\ : 0,2 : az equal to 0.16 : 0.36 : 0.48 
and a phase difference of 180° between the 100 Hz component 
and the two other components. 

Using 2:3:4 driving we indeed observed a perfect 8- 
fold quasi-pattern. In Figure |43j we show images and 
power spectra of this state at two temporal phases. This 
state is subharmonic in time and can be observed in the 
region where 10-fold quasi-patterns and 8-fold distorted 
quasi-patterns were observed for the 2/3 experiments de- 
scribed above. Interestingly, this state was exceedingly 
stable and existed within a single domain over a wide 
range of parameters (note the sharp peaks in Figure O) . 
This is in sharp contrast to the distorted 8-fold state 
shown in Figure [42 which existed in both a narrow range 
of parameters and, as evident in its diffuse spectrum, had 
a tendency to break up into domains. It is possible that 
the third frequency allows the 2kR mechanism to form 
the quasi-pattern through the relaxation of the ratio be- 
tween k\ and k 2 and thereby the angle between them. 
This state was observed in a 50/75/100 Hz experiment 
with ratios of driving amplitudes a\ : a 2 : 03 equal to 
0.16 : 0.36 : 0.48, where the single-frequency critical ac- 
celerations for these frequencies are, respectively, 2.56g, 
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4.40g and 6.91g (yielding ratios of a 0.18 : 0.32 : 0.49) 
and phase differences <j)\ = 0°, <j>\ = 180°. Although the 
third frequency acceleration is not small here, it is still 
below the critical value for single-frequency excitation 
at 100Hz and wavevectors corresponding to the 100 Hz 
component were not observed. 



E. Transition states to rhomboids 







Figure H shows an expanded view of the 2kR phase 
space around \ c (see Figure H (right)). The region shown, 
lies in the subharmonic region of phase space (x > Xc) 
where square patterns form at threshold. For \ close to 
Xc the transition region and the 2kR region is relatively 
narrow compared to the transition region for larger values 
of X- 
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Driving amplitude (g) 50 Hz. 
FIG. 44. An expanded view of the phase space described 
in Figure H (right). Transition states between square and 2kR 
states are formed by a superposition of a square state with 
one (two) sets of 2kR triad vectors. Between the 2kR and 
transition states, a modulated 2kR state exists in which an 
additional fc2 component parallel to the k\ component exists. 

The transition from the square symmetry region to the 
2kR state is perhaps more interesting than the transition 
between squares to square 2MS discussed in (subsection 
VIC). While the square-2MS transition occurs through 
the formation of an "additional" set of wavevectors (of 
different wavenumber) that combine with each of the pri- 
mary wavevectors that initially formed the squares, the 
square-2fci? transition has a qualitatively different na- 
ture. Here, the basic square symmetry is not only bro- 
ken, but is actually replaced by a pattern of completely 
different symmetry. In 2kR states, one of the two fc 2 
wavevectors that are initially perpendicular to each other 
in the square state is replaced by a fc 2 wavevector whose 
orientation forms the angle 8 that is determined by the 
magnitude of the k\ wavevector which defines the 2kR 
rhomboidal pattern. 



FIG. 45. A typical time sequence of the transition state is 
shown with the corresponding power spectra (bottom). The 
wavevectors (k2 a , A?2&) generating the squares appear together 
with a third wavevector, ki c that is oriented at an angle of 
6 — 41° to k-2b- The above experiment was made for system 
parameters of lu /{2tt) - 27Hz, v = 22.cS, (p = 0°, x = 69° 
and h = 0.2cm. 

A typical "transition" state at different temporal 
phases is shown in Figure Ea. The transition occurs 
through the formation of one or two additional wavevec- 
tors of magnitude k%, (&2c in Figure Hqa). These new 
wavevectors are aligned at the 2kR resonant angle 6 with 
respect to the fc 2 wavevectors (fc 2a arL d &2&) which form 
the squares. Although additional fc 2 vectors could, in the- 
ory, form 2kR triads with all of the original fc 2 wavevec- 
tors, empirically, we find that only one of the initial fc 2 
directions is selected. This is possibly due to the fact 
that the self-interaction of the harmonic k\ wavevectors 
prefers a hexagonal (rather than a square) arrangement. 
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FIG. 46. The transition states between square and 2kR 
symmetry do not usually form a global state, (a) and (b) show 
different phases of a typical transition state that is formed by 
three domains. In (a) white circles note areas where two 2kR 
domains intersect. The power spectrum of (a) is shown in (c) 
and the wavevectors are noted in (d) by (squares) the primary 
square wavevectors, (circles) 2fc7?-resonant ki wavevectors, 
and (diamonds) the k\ wavevectors. Since the power spectra 
reflect the entire pattern, three separate sets of 2kR triads can 
be observed in the same spectrum. The system parameters 
here are lu /(2tt) = SOHz, v = 23cS, h = 0.2cm, x = 71.4° 
and<£ = 0°. 

Thus the transition pattern that is formed is a superpo- 
sition of a k% square state with a 2kR state. All other 
wavevectors seen are formed by secondary interactions of 
these wavevectors. Generally, the transition state breaks 
the symmetry of the square pattern in a single direc- 
tion. Since the transition state is usually not global, two 
domains can form and at their common border one can 
observe structures that retain the four-fold symmetry, as 
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indicated by the white circles in Figure |4qa. 

In the near vicinity of Xci the transition patterns ap- 
pear as in Figure [45], and their rhomboidal character is 
apparent. Farther away from Xc the appearance of the 
transition state is more similar to squares, as the orthog- 
onal fa wavevectors are more dominant. 
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FIG. 47. "Modulated 2kR" states (left) are variants of the 
2kR states formed when (see the power spectrum (right)) an 
additional wavevector of magnitude fa, fab, parallel to fa is 
generated. The system parameters are u>o/{2n) — 22 Hz, 
v = 23 cS, <f> = 0° and h = 0.2 cm. 

For x > 70° we observe a variant of the 2kR state, 
shown in Figure |47J, which forms between the "transi- 
tion" and the 2kR regions (see Figure |44|). This state, 
which we call a "modulated" 2kR, consists of a modula- 
tion of the regular 2kR triad by an additional wavevec- 
tor of the magnitude and temporal behavior of the larger 
frequency component. This state is formed by the super- 
position of a 2kR state with an additional wave vector of 
magnitude fa that is oriented parallel to the k\ wave vec- 
tor (see Figure [47] (right)). This wavevector retains the 
time dependence of the fa mode and generates a spa- 
tial and temporal modulation, similar to that induced 
in the SSS mode, in the fa direction (x in Figure 34 1). 
A possible mechanism that can generate a "modulated" 
2kR state is a 2MS resonance between the two collinear 
wavevectors fa and fa that produces a fa wavevector 
along the same direction. This is verified by comparison 
with the square 2MS pattern obtained for the same pa- 
rameters with <f> = 90° (see Figure B8e). This mechanism 
is supported by our observations of the "modulated" 2kR 
states for only relatively low values of ojq for which the 
wavevectors fa, fa and fa, forming the 2MS state, are 
nearly parallel. 



VIII. ADDITIONAL RESONANT STATES 

A. States observed for odd/odd parity driving 

The work of Silber et al. |53| indicates that the pos- 
sible three-wave resonant interactions between excited 
and damped modes depend on the parity of the driv- 
ing. For odd/odd driving no 3-wave resonant interac- 
tions between fa and fa are expected since both modes 



are temporally subharmonic (see Sec. ID 1). When con- 
sidering four wave interactions, however, these restric- 
tions are no longer valid. Below, we describe two ex- 
periments conducted with odd/odd driving that indeed 
show that 4-wave resonant interactions can be found 
both in the vicinity of Xc where one mode is strongly 
excited and the other only weakly damped, and far from 
Xc where the interactions involve an excited mode and a 
single strongly damped mode corresponding to the sub- 
harmonic frequency, loq/2. 
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FIG. 48. A 5/7 ratio yields both the 2MS resonance (a) 
and a four-wave resonance (b). In the 2MS power spectrum 
(a right) the three vectors are nearly parallel obeying the res- 
onance fa = fa — fa . In (b right) the power spectrum reveals 
a qualitatively different resonance: fa + fa = fa + k' 2 = fa. 
This resonance produces additional vectors such as fa that 
are sums and differences of the original wavevectors involved. 
The experimental parameters are 55/77 Hz driving frequen- 
cies, v = 23 cS, and h = 0.2 cm. In (a) 4> = 90°, \ = 72.8° 
and in (b) = 0°, x = ™.8°. 

For x > Xc a 5/7 frequency ratio can produce a square 
2MS state in the vicinity Xc- This, however, is not the 
only resonantly locked pattern that can be observed at 
this driving ratio. Mixing angles typically 3-5 degrees 
beyond Xc results in the formation of the state shown 
in Figure J48p. The spatial spectrum of this state indi- 
cates that it is formed by a qualitatively different mech- 
anism. The original square pattern, formed by two or- 
thogonal fa wavevectors, is broken by an additional pair 
of fa wavevectors, whose orientation is determined by the 
four-wave resonance condition: fa + k[ = fa + k 2 - This 
resonance produces additional vectors that are sums and 
differences of the original wavevectors involved. States 
similar to these have been previously observed in ver- 
tically oscillating convection experiments [ p2[ , where 3- 
wave interactions arc forbidden. 
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FIG. 49. In an experiment with 3/5 (60/100 Hz) driving, 
patterns with square symmetry (or slightly rhombic) domi- 
nate for mixing angles x < Xc- The primary wavenumber is 
of magnitude fci (circles in a,b,c). When the driving is in- 
creased, the primary symmetry is broken by a small wavevec- 
tor of magnitude k^ (enclosed in triangles). These states are 
observed for parameters v = 23 cS, and h = 0.2 cm both 
at mixing angle close to \c ( a ) X — 69.6° and far from Xc 
(b) x — 48.1° and (c) x — 58.7°. A possible mechanism for 
the large scale symmetry breaking is a resonant four-wave 
interaction. The power spectrum (a bottom) reveals the ex- 
istence of additional vectors of magnitude fci (enclosed by 
squares). The interaction of k\ c with the sum of the origi- 
nal vectors k\ a + kib (enclosed by diamonds) results in a new 
vector ks — k\ c — (k\ a + fcib) whose scale is consistent with 

In odd/odd driving both sides of phase space exhibit 
subharmonic temporal response. Unlike the 5/7 driving 
that produces hexagons for x < Xc (despite the subhar- 
monic temporal response), experiments performed using 
3/5 driving with luq = 20 Hz yield square symmetric 
patterns for x > Xc an d nearly square patterns (rhom- 
boids with an angle of 84°) for x < Xc- Here, a state 
exhibiting a different mode of 4-wave coupling far from 
Xc is presented in Figure |4Ej. Like the SSS states, these 
patterns result from a symmetry breaking bifurcation for 
Xc + 20° > x > Xc- The power spectrum reveals that 
in addition to the dominant pair of orthogonal wavevec- 
tors, (fci a and ku, in Figure |49[ ) an additional wavevector 
of the same magnitude, fci c , is created along the bisector 
of k\ a and ku,. The vector sum of these three wavevec- 
tors of magnitude k\ produces a smaller wavevector of 
magnitude k^ = k\ c — (k\ a + ku,). The scale of k$ is 
consistent with k(ujo/2), as determined by the linear dis- 
persion relation. Thus, as in the case of 2MS states, a 
symmetry-breaking slaved mode, k% is excited by a non- 
linear resonance. 



B. States satisfying more than one resonance 
condition 



We have observed a number of cases where states which 
satisfy more than a single resonance condition were se- 



lected by the system. These states are generally stable in 
a relatively wide range of phase space. Here we present a 
number of examples of such multiply-resonant nonlinear 
states. 

We have seen that SSS-I states result from a primary 
hexagonal symmetry broken by a wavevector of size k c /2. 
A similar mechanism can occur for square or rhomboid 
patterns. In Figure |50] we show a spatially subharmonic 
state, obtained using 5/8 driving, where a rhomboid, is 
formed by two wavevectors of magnitude k\ (fci and k[ in 
Figure |5(j) via the additive 2kR resonance: ki + k[ = k 2 - 
The spatial period of the rhomboid is doubled by the 
appearance of two new wavevectors of size ki/2, (k\j2 
and k[/2 in Figure p0). This state was observed for a 
mixing angle x = Xc— 1°< An additional spatial resonance 
is apparent. The vectors, k\ — k[ + ki/2 and k\ — k[ — k[/2 
both form vectors of length ki (see k 2a and kib in Figure 
|50| ) . This "extra" 4-wave resonance may be the result of 
spatial mode-locking. It may be possible that this extra 
resonance causes the selection and resultant stability of 
this state. 
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FIG. 50. An image (left) and power spectrum (right) of 
a spatially subharmonic state where a rhomboid is formed 
by ki + k[ — ki where k\ (fe) are generated by the 5u>o 
(8a;o) frequencies. The pure rhomboid symmetry is broken 
by two spatially subharmonic vectors, fci/2 and fci/2. The 
vectors kia and kib are formed by additional resonances: 
kia — k\ — k[ + fci/2 and kib — k\ — k[ — k[/2. This state 
appears for 50/80 Hz driving for v — 23 cS and h = 0.2 cm 
for X~Xc- 1°- 

Another example of a state satisfying two resonance 
conditions is presented in Figure KU This state, which 
is only observed in 1/2 driving experiments, is both a 
hexagonal 2MS state as well as an SSS-II state. The driv- 
ing ratio of 1/2 is unique since U3 (given by u>3 = L02 — u>i) 
is equal to ui\. Thus, the wavenumber &3, excited by 0*3, 
coincides with the wavenumber k\ , excited by uj\ . There- 
fore n] , a resonant triad is formed involving only the two 
critical wavenumbers k\ and ki- In this case SSS-II and 
2MS states coincide for 20 Hz< c^o <50 Hz. The reso- 
nance that is formed is identical to the SSS-II resonance 
found for odd/even driving in the harmonic region (see 
Sec. IV B), where the wavevectors of magnitude q and 
K of the SSS-II state are replaced by wavevectors of re- 
spective magnitude k\ and k^. The resonance condition 
for the 1/2 experiment can be written as fci + k[ = k-z. 
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The SSS-II type resonance is possible since, for the sys- 
tem parameters used, the ratio fo/fci is close to y/3. The 
same resonance was also observed by Wagner et al. for a 
two mode subharmonic-harmonic interaction in a single- 
frequency, lubrication-limit experiment E2]. In their 
experiment, for special values of h and w, two modes 
with temporal responses of lu/2 and u> and with spatial 
wavenumbers k\ and ki-, are excited by single-frequency 
forcing. 
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FIG. 52. A state formed by simultaneous multiple reso- 
nances obtained with 2/3 driving with the same parameters 
used to form a 2kR state at <f) = 0, cjo/(2tt) = 20 Hz, v — 23 
cS, x — 59°, and h — 0.2 cm, but with a phase of <f> — 117°. 
The wavevectors of magnitude k$ are formed both from the 
resonance fo — k' 2 — fc 3 ' as well as the resonance k?, — k\ — hj, . 
In addition, the condition satisfying the DHS structure is sat- 
isfied and the hexagonal sublattice formed by the k$ wavevec- 
tors spans the entire lattice, ki, fe, and kz are, respectively, 
noted by squares, circles, and diamond symbols. 






FIG. 51. Images (a-c top) and the corresponding power 
spectra (a-c bottom) of different temporal phases of a SSS-II 
(2MS) state observed for 1/2 driving in the harmonic region 
of phase space at system parameters of v — 23 cS, 40/80 Hz 
and h = 0.155 cm. The pattern was observed in the vicinity 
of Xc for x > Xc- This state resembles SSS-II states observed 
for 2/3 driving although here a typical super-hexagon cell is 
it/3 symmetric whereas the SSS-II state generally has n/6 
rotational symmetry. 

In our last example, we present an example of a state 
where multiple resonance conditions are simultaneously 
satisfied. When this occurs, we find that the stability of 
a pattern is significantly enhanced. This state contains 
features of most of the states described above. This single 
state includes a 2MS resonance, a DHS resonance (with a 
1/V7 sublattice), an SSS-II type resonance and an "ad- 
ditive" and "subtractive" 2kR resonance. This state was 
observed for 2/3 driving (o>o = 20 Hz) with = 117° and 
is shown in Figure |52J. 

As can be seen in the power spectrum (Figure |32j right) 
a double lattice with wavenumber k^ at an angle of 22.2° 
produces the familiar DHS structure described in Section 



[II. However, in this case the smallest inner hexagon 
has a magnitude of k^, which is associated with w 3 = 
LU2 — u>i frequency, while the second inner hexagonal set 
of wavevectors has a magnitude of k\ . This resonance is 
very stable and exists in a single domain. This stability 
is perhaps due its multiply-resonant nature. In Figure 
p3| a time sequence is shown of the different phases of this 
state is presented. 




FIG. 53. Different temporal phases of the state described in 
Figure M taken for constant values of the driving parameters. 
Circles of radii k\ (middle circle), k-z (outer circle) and £3 
(inner circle) are drawn in (b). At different temporal phases, 
different scales appear dominant. 



IX. CONCLUSIONS 

The work described in this paper provides, for the first 
time, a partial, but coherent, experimental picture of 
both the nonlinear states generated by two-frequency 
forcing as well as their domains of existence in phase 
space and the nonlinear interactions that generate them. 
As demonstrated above, the space of nonlinear patterns 
formed by two interacting unstable modes is very rich. 
Our understanding of the types of structures and their se- 
lection is just beginning. Predicted nonlinear three-wave 
resonances (3(J, however, appear to govern nonlinear pat- 
tern selection only for the simplest ratios (2/3 and 1/2). 
Recent theoretical advances made by Silber et al. in the 
analysis of two-frequency Faraday system in the vicin- 
ity of the co-dimension 2 point |30|,[33| suggest that at 
least some of the features discovered by our experiments 
can be reproduced by amplitude equations derived using 
the quasi-potential approximation (e.g. pT{|). In addi- 
tion, recent work |2^,^,^1[ indicates that many of the 
superlattice states found far from Xc may be qualitatively 
understood as representations of invariant subgroups of 
broken hexagonal symmetry. A quantitative theory that 
describes the important regime of two or more concur- 
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rently unstable modes for arbitrary height and viscosity 
has still to be developed. 

Some of the major results of this work are summarized 
below: 

• The temporal symmetry of the driving determines 
which resonant interactions can occur between the 
primary excited modes. A number of distinct three 
and four wave interactions between excited modes 
and between slaved and excited modes were exper- 
imentally observed. 

• Factors such as the dimensionless dissipation in the 
system and the driving phase <fi can play an impor- 
tant role in the nonlinear interaction mechanism 
selected by the system. (f> is a convenient param- 
eter for comparison with theory, as, in contrast to 
the system's dissipation, its value does not play a 
role in the validity of the theoretical (e.g. IO]) ap- 
proximation used. 

• Symmetry-breaking can often occur via modes that 
nonlinearly couple to the original set of modes. The 
symmetry-breaking modes tend to be arranged 
in invariant subgroups of the original symmetry 
group. 

• In many cases, slaved modes can be non-linearly 
excited by the parametrically amplified linear 
modes. The availability of slaved modes (deter- 
mined by the ratio of the driving frequencies) is im- 
portant in the selection of the final nonlinear states. 
This provides a nontrivial selection mechanism for 
nonlinear states. 



• A theoretically predicted mechanism [ (25|]26| for 
producing quasi-patterns was for the first time ex- 
perimentally observed. By tuning the system pa- 
rameters to satisfy a resonance condition 2n-fold 
quasi patterns can be produced for any desired n. 



We believe that both the states and nonlinear mech- 
anisms described in this work should be of general im- 
portance to a wide class of parametrically-driven nonlin- 
ear systems. Such systems include parametrically-driven 
fluid systems, non-linear optical systems [f^|64|]5l| ,65,44|, 
nonlinear wave interactions in superfluid helium [56,37|, 
magnetically-driven ferrofluids |45[ | , and possibly nonlin- 
early coupled mechanical systems. Although the behav- 
ior of systems driven by 2-frequency forcing is, itself, 
important. Understanding the spatio-temporal behavior 
of such systems is but a first step in understanding the 
general behavior in space and time of nonlinear systems 
driven by multiple frequencies. The work presented here 
is an important building block on the road to understand- 
ing these more complex systems. 

We gratefully acknowledge the support of the Israel 
Academy of Sciences (Grant no. 203/99). 
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